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Abstract: We discuss various aspects of resummed chiral perturbation theory, which was devel- 
oped recently in order to consistently include the possibility of large vacuum fluctuations of the 
ss-pairs and the scenario with smaller value of the qq condensate for Nf = 3. The subtleties of 
this approach are illustrated using a concrete example of observables connected with ttt; scattering. 
This process seems to be a suitable theoretical laboratory for this purpose due to its sensitivity to 
the values of the O(p^) LEC's, namely to the values of the fluctuation parameters L4 and Lq. We 
discuss several issues in detail, namely the choice of "good" observables and properties of their bare 
expansions, the "safe" reparametrization in terms of physical observables, the implementation of 
exact perturbative unitarity and exact renormalization scale independence, the role of higher order 
remainders and their estimates. We make a detailed comparison with standard chiral perturbation 
theory and use generalized well as resonance chiral theory to estimate the higher order 

remainders. 
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1. Introduction 



As it is well known, at the energy scales E <^ Ah ^ IGcV the physics of QCD is nonperturbative 
and governed by chiral symmetry (xS") SU{Nf)LxSU{Nf)ji. This global symmetry is present on 
the classical level within the QCD with Nf massless quarks (in the chiral limit of QCD) and on 
the quantum level there exist strong theoretical (for Nf > 3) and phcnomcnological arguments for 
spontaneous symmetry breakdown (SSB) of according to the pattern SU{Nf)LxSU{Nf)ii 
SU{Nf)v- Due to the confinement, quark and gluon fields do not represent appropriate low energy 
degrees of freedom within the above mentioned energy range; the relevant degrees of freedom 
correspond to the lightest colourless hadrons in the QCD spectrum. As far as the Green functions 
of quark currents are concerned, it is possible to obtain a general solution of the chiral Ward 
identities in terms of the low energy expansion. This expansion can be organized most efficiently 
using the methods of effective field theory corresponding to the low-energy limit of QCD with Nf 
light quark flavours which is known as chiral perturbation theory (xPT) 0, ^ ^ . xPT describes the 
low energy QCD dynamics in terms of the lightest (A^| — l)-plet of the pseudoscalar mesons identified 
with the Goldstone bosons (GB) of the spontaneously broken chiral symmetry which appear in the 
particle spectrum of the theory as a consequence of the Goldstone theorem. In the chiral limit 
these pseudoscalars are massless and dominate the low energy dynamics of QCD. They interact 
weekly at low energies E <C Ajy, where ~ IGeV is the hadronic scale corresponding to the 
masses of the lightest nongoldstone hadrons. This feature of the GB dynamics enables systematic 
perturbative treatment with the expansion parameter (E/Ah)- Within the real QCD the quark 
mass term jCj ^^^^ breaks explicitly and the Goldstone bosons become pseudogoldstone bosons 
(PGB) with nonzero masses. Though nif ^ 0, for ruf <^ Ah the mass term C^j'^^g can be 
treated as a perturbation. As a consequence, PGB correspond to the lightest hadrons in the QCD 
spectrum^ (identified with 7r°,7r^ for Nf = 2 and , , , ,?] for Nf — 3) and the 

interaction of PGB at the energy scale E ^ Ah continues to be weak. Because Mp < Ah, the 
QCD dynamics a,t E ^ Ah is still dominated by these particles and the effective theory provides 
us with a simultaneous expansion in powers of {E/Ah) and (m//Aij). The Lagrangian of xPP 
can be constructed on the basis of symmetry arguments only; the unknown information about the 
nonperturbative properties of QCD are hidden in the parameters known as low energy constants 
(LEC)I^, 1^. These are related to the (generally nonlocal) order parameters of the SSB of x^, the 
most prominent of them are the Goldstone boson decay constant Fq and the chiral condensate^ 
Bo = J:/F§ where E = -{uu)o . 

To be more precise, 7V/-flavour xPT is in fact an expansion in rrii, around the SU{Nf)L x 
SU{Nf)ii chiral limit rrii = 0, i < Nf, while keeping all the other quark masses for i > Nf at their 
physical values. Because niu.d are much smaller not only in comparison with the hadronic scale 

^The PGB masses Mp can be expanded in the powers (and logarithms) of the quark masses starting from the 
linear term and therefore vanish in the chiral limit. 

^The parameter _Fo is however more fundamental in the sense that Fq is both necessary and sufficient condition 
for SSB, while {qJqf)o 7^ corresponds to the sufficient condition only. (The lower index zero means here the chiral 
limit.) 
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A//, but also in comparison with the intrinsic QCD scale Aqcd^ the two-flavour xP^I" is expected to 
produce well-behaved expansion corresponding to small corrections to the SU{2)l x SU{2)u chiral 
limit. 

The strange quark mass on the other hand, though still small enough with respect to Ah to 
be treated as an expansion parameter within the three-flavour x^^^ (relating real QCD with its 
SU{3)l X SU{3)n chiral limit), is of comparable size with respect to Aqcd- This fact, besides the 
expected worse convergence of the three-flavour xPT, might also have interesting consequences for 
the possible difference between the Nj — 2 and Nf = 3 chiral dynamics. As discussed intensively 
in a series of papers IH ^, ||, 0, ||] , rn^ < Aqcd suggest, that the loop effects of the vacuum ss 
pairs are not suppressed as strongly as it is for the heavy quarks and might enhance the magnitude 
of the Nf — 2 chiral order parameters relatively to their Nf = 3 chiral limits. This applies mainly 
to Fo{Nf) and Y.{Nf) — FQ{Nf)Bo{Nf), which should satisfy paramagnetic inequalities |Q| 

S(2) > E(3) = lim I](2) 

^^o(2) > Fo(3) = lim Fo(2). (1.1) 

rUs — ^0 

The leading order difference between the two-flavour and three-flavour values is proportional to rus , 
with coefficients measuring the violation of the OZI rule in the 0^+ channel, e.g. 

i:{2) = + m,zl + . . . (1.2) 



(see B for details) where 



+0 



Zi = lim / (rx{uu{x)ss{0))c (1.3) 



and analogously for Fq. The fluctuation parameter Zl is related to the LEG Lg{^) (and L\{ii) 
for Fo) of the three-flavour xPT • As discussed in [|[ g §, these parameters might be 
larger than their estimate based on the large Nc expansion, provided Nf = 3 is close to the critical 
number of light quark flavours Nj"'^, for which the chiral symmetry is restored. Available estimates 
vary widely, some indicate a larger number JV^"* - 10 - 12 for A^^ = 3 (lO[ |ll|, |l|, while other 
approaches |l^, 14 1 and lattice calculations jl5[ |6[ ^ |l8[ |l^ discuss a possibly much lower value 



Nf"* < 6. Provided the scenario of large vacuum fluctuations takes place, the second term in (1.2) 
(called the induced condensate in ^ |2^ ) can be numerically comparable with the flrst term and 
the three-flavour condensate S(3) could be substantially smaller than the two-flavour one, the value 
of which is experimentally accessible in the recent experiments. Analogous reasonings apply to the 
relationship of -Fo(2) and -Fo(3). 

These effects could possibly have strong consequences for the organization of the chiral expan- 
sion in the Nf = 3 case [§, §|, @, D . Let us remind that the general form of the Lagrangian of xPT 
is 

/: = ^ (1.4) 

m,n 

where 

k 

with leg's (7^™'"-' and independent set of the operators O^™'"-* = 0(5"*™^). 

In order to be able to treat the double expansion consistently, it is necessary to assign a single 
integer parameter called chiral order to each term = 0{d"^m^) of the effective Lagrangian. 

The terms Ck with chiral order k are then called 0{p^) terms. Obviously, d — 0{p). The matter 
of discussion might be, however, the question concerning the chiral power of ruf. This question is 
intimately connected to the scenario according which the SSB of xS is realized. 
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The standard scenario Ig, ^ corresponds to the assumption, that the SSB order parameters 
E(7V/) and Fo^Nf) are large in the sense, that the ratios 



(where and m — (m„ + md)/2) and 



^ -TT 



are close to one. Because = 0{p'^), it is then natural to take ruf ^ 0{p^), i.e. k = m + 2n. This 
results in the standard (SxPT in what follows). This scenario seems to be experimentally 

confirmed ||2l[ for Nf — 2; the recent analysis of the data yields |^ 

X(2) = 0.81 ±0.07, Z(2) = 0.89 ±0.03. (1.8) 

The 0(p2) Lagrangian |, | 

^2^^ {{df,U+d''U) + 2Bo{U+M + M+U)) (1.9) 

gives E(2)lo = S(3) = BqFq at the leading order, thus postponing the difference S(2)— S(3) to 
higher orders. The same is true for the parameters Fg^Nf). Let us also note that the quark mass 
ratio r — mg/rh is not a free parameter here'^, at the leading order one has 

(1,10) 

TT 

An alternative way of chiral power counting for Nf = 3 is the generalized (GxPT) 
|2^, 11^, |2^, |2^, , originally designed to treat the scenario with small quark condensate X{3) 
and to take the quark mass ratio r as a free parameter. In the case X(3) ^ 1 it is natural to take 
TO/ = 0{p) and Bo — 0{p), this means k — m + n. In contrast to SxPT, there are also odd chiral 
orders and the 0{p^) Lagrangian contains additional terms which are O(p^) within the standard 
chiral counting^ (see e.g. |2^, |2^): 

>C2 = ^ {{d^,U+d^U) + 2Bo{U+M + M + U) + Aa{{U+Mf + {M + Uf) 

+ Z^{U+M - M+Uf + Z^{U+M + M+Uf) . (1.11) 

For the condensate S = ~{uu) we get at the leading order for Nf — 3 

Slo BoF^ + Z^{2m + to,) = E(3) + Zo^(2to + to,) (1.12) 

and therefore 

I]lo(2) = S(3) + Zo^to,. (1.13) 

This allows the difference S(2)— 5](3) to appear already at the leading order, consistently with the 
small S(3) scenario. The next-to-leading order Lagrangian 0{p^) 

£3 = ^ {^^{df.U+d^'UU+M + M+U) + i,{d^U+d''U){U+M + M+U) + . . . j (1.14) 



''On the contrary, the value of r is usually taken as an input in standard 0{p^) fits, see e.g. |23| and references 
therein. 

■^Effectively the generalized chiral power-counting means partial resummation of these terms. 
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(where the eUipses stand for the additional terms which are of the order 0{p^) in SxPT) gives rise 
to the Nf = 3 relation 

FInlo = ^^0 (3)(1 + 2e(m, + 2m) + 2mO, (1-15) 
which implies that the difference Fq{2) — -Fq (3) is treated as an effect of the next-to- leading order 

F2(2) = F2(3)(f + 2fm,). (1.16) 

Therefore, neither SxPT nor GxPT can accumulate the case of large fluctuation parameter ^ and 
the ratio Z{3) ^ 1 at the leading order. 

Quite recently, a consistent method of handling the case X{3), Z(3) < 1 was proposed [g |, @, 
H, |9| . Instead of changing the chiral power counting, it is based on a more careful manipulations with 
the chiral expansion. As it was discussed in the above references, the case ^(3), Z{3) ^ 1 could 
significantly influence the properties of the chiral expansion inducing instabilities of the perturbative 
series corresponding to the observables, which cannot be linearly related to the QCD correlators 
(such as the ratios like PGB masses, scattering amplitudes etc.). For such quantities, one should 
not perform a perturbative chiral expansion of the denominators but rather keep the ratios in 
a nonperturbative "resummed" form. The possibly large vacuum ss pair fluctuations are then 
parameterized in terms of X{3), Z{3) and r and treated as free parameters. We return to the 
detailed formulation of this recipe in the next section. 

The aim of this paper is to illustrate the "resummed" form of the chiral expansion with spe- 
cial attention to its formal properties and to the details and subtleties of the general procedure. 
Motivated by our preliminary results on nrj scattering within the GxPT we have chosen the 
observables connected with this process as a concrete example which seems to be sensitive to the 
deviations from the standard assumption ^(3), Z(3) ^ 1 (note that some recent phenomenological 
studies suggest a possibility of ^(3) ~ 0.5, cf. ||3^, ^ and Q). Also, from the phenomeno- 
logical point of view, the off-shell TTrjTrr]* vertex is a necessary building block for the non-resonant 
part of the amplitude for the rare decay rj — > 7r'^7r'^77. Preliminary estimates within GxPT |3^, |3^ 
suggest, that the effect of deviation of this off-shell vertex from the standard case might be at least 
in principle observed. The details will be presented elsewhere fs^ . 

The amplitude of tti] scattering was already calculated within SxPT to 0{p^) (and within 
the extended SxPT with explicit resonance fields) in the paper |Q, where the authors presented 
prediction for the scattering lengths and phase shifts of the S, P and D partial waves. We quote 
here their 0{p'^) results for the S- and P-wave scattering lengths (in the units of the pion Compton 
wavelength): Qq^^^ = 7.2 x IQ-^ and af^^^ = -5.2 x 10"^. 

The paper is organized as follows. In Section 2 we recapitulate the motivation for the resummed 
version of xPT and the construction of the bare expansion of "good" observables. We make a 
detailed general discussion, connected with the four-mcson amplitude, of the strict chiral expansion, 
the dispersive representation and matching of both approaches with stress to the reconciling of 
exact perturbative unitarity and exact renormalization scale independence in Section 3. Section 4 
is devoted to the general properties of the nrj scattering amplitude. We discuss the kinematics, the 
definition of suitable "good" observables, the dispersion representation of the amplitude and the 
construction of the bare expansion. Various possibilities of its reparametrization are described in a 
detailed way in Section 5. The numerical illustration of the particular variants is made in Section 
6, where we also numerically illustrate the subtleties of the construction of the bare expansion. 
We recapitulate the results of the standard variant of xPT and compare them with the resummed 
approach. We concentrate on the dependence on the LEC's as well as on the sensitivity to the 
higher order reminders and make an attempt to estimate their values using a matching with GxPT 
and a simple version of resonance chiral theory. Section Q contains the summary and conclusions. 
Some technical details are postponed to the appendices. 
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2. Resummation of the vacuum fluctuations - motivation and basic nota- 
tion 



As we have mentioned in the Introduction, the potentiaUy large vacuum fluctuations of the ss pairs 
might result in the instabilities of the chiral expansion, which originate in the possibility that for 
some observables the next-to-leading order correction could be numerically comparable with the 
leading order one. As discussed in this could generally cause problems with the convergence of 

the formal chiral expansion. Nevertheless, at least for some carefully defined "good" observables, it 
is natural to presume some sort of satisfactory convergence properties. Such "good" observables are 
assumed to be those which can be obtained directly from the low-energy correlation functions in the 
domain of their analyticity far away from singularities and which are related to the corresponding 
correlator linearly ^J. Typical examples are the squares of the PGB decay constants Fp , the 
products FpMp where Mp are the PGB masses and also the subthreshold parameters which can 
be derived from the products A Yii^i where A is the PGB scattering amplitude 1-1-2^3 + 4. 
Let us write the expansion of such a "good" observable G in the form of a (carefully defined) hare 
expansion |8j as^ 

where G^^^ = g^'^\Fo, Bo,mq) and G'^^ = q'^^^Fq, Bq, niq, Li, Mp) correspond to the sum of the 
leading and next-to-leading order terms respectively and the renormalization scale independent 
quantity Sq accommodates the higher order remainders. 

As a terminological note, in what follows we use the term strict chiral expansion for an un- 
modified expansion in terms of the LEC's strictly respecting the chiral orders. The bare expansion, 
though still expressed in terms of LEC's, accumulates some modifications dictated by physical 
requirements. It is the bare expansion which is assumed to be globally convergent. 

For a "good" observable it is then assumed 

\Sg\ « 1 (2.2) 



as a natural assumption. This property of the bare chiral expansion (2.1 ) is called global convergence 
in 1^, p. Note however, that the validity of the inequality ( |2.2| ) might depend on the definition of 
the reminder Sq which is not fixed unambiguously and might differ according to the calculation 
scheme in use. We will comment on this point later on. 



The above mentioned possible instability in (2.1) appears when G(2) G'^'^\ i.e. ATg $ 1, 
where 

G(2) 

^G = — . (2.3) 

Such an instability manifests itself in the expansion of the observables depending on G nonlinearly 
For instance, for a ratio of two "good" observables G and G formally expanded in the form 



G /G(2)\ /G(2)\/gW g'W\ G 



W~{G^J^yG^J \g(^'g^ j^a^''/'''' ^^-^^ 

we get for the remainder Sq jQi 

^G/G' - + ^ - (2-5) 

For Xq' < 1 this might be numerically large even if both \6g\, \Sq' \ were reasonably small. In this 
sense, a ratio of two globally convergent observables need not to be necessarily globally convergent 



^Here we tacitly assume the standard chiral power counting. Analogous expansion could be written also for the 
generalized case. 
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too. It should be therefore much safer not to expand such "dangerous" observables and rather write 
the ratio in the "resummed" form 

G _ Gt^) + G'^) Gj 

G " G'(2) + G'W + G^'^G/G' ■ (^-(5) 
The relation ( |2.6| ) is an exact algebraic identity provided we keep explicitly the remainder 

In this case (5^^^' remains for |(5g|j I<5g'I ^ ^ under numerical control. 

Of course, only the fact that the bare expansion of some observable is not globally convergent 
does not necessarily correspond to the collapse of the convergence, because the next-to-next-to- 
leading order G^^^ can saturate the series in such a way that the next-to-next-to-leading remainder 

Q^NNLO q(2) _ ^(4) _ g(6) (2.8) 

is reasonably small. Namely this is the usual assumption behind the O(p^) calculations. Violation 
of the global convergence property here means merely that the 0{p^) contribution have unnatural 
size, i.e. G'^-* < G'^' -f- G^'^\ This could, however, destabilize the 0{p^) chiral expansion of ratios 
in the way similar to that discussed above. 

Provided we allow the expansion of the "good" observables only, we are also pressed to modify 
the next step leading from the bare expansion to the usual output of the x^T, consisting of 
a reparametrization of the expansion by expressing some of the LEC's in terms of the physical 
observables such as masses and PGB decay constants. This step converts the series into an expansion 
in powers and logs of the (squared) PGB masses instead of quark masses. To achieve this, it is either 
necessary to invert a bare chiral expansion of some observable (in the case of the 0{p^) LEC's) or 
to use an observable which might be generally a "dangerous" one. Let us briefly discuss the first 
case. Schematically, suppose that some O(p^) LEG Go {e.g. Fg) just corresponds to the leading 
term G^^) of the expansion of the observable G. Then we can write an algebraic identity 

Go = G - G^-^) (Go) - G6g, (2.9) 

where we explicitly point out the dependence of the next-to-leading term on Go. To convert this 
expansion and express Go by means of the series in G one substitutes G for Go on the right hand 
side. This defines a new remainder 

Go=G-G(4)(G)+Go<5go, (2.10) 

for which we get 

A ^-Xg , 1 GW(G) 

dGo = y + ^ T^, • (2-11) 

Ag Ag Cj 

This could cause an instability of the converted expansion for Go in terms of G for Xg ^ 1 even 
if the relative size of the next-to-leading order G^^\G)/G is reasonably small, irrespective of the 
condition for global convergence \5g\ <C 1. 

On the other hand, suppose that some 0{p'^) constant Gi coincides with the next-to-lcading 
term G^'^K In this case we have an algebraic identity for Gi 

Gi = G - G(2) - G(5g (2.12) 

and the remainder here is perfectly under control, provided G has a globally convergent bare 
expansion and we do not re-express G*-^' in terms of physical observables (i.e. provided we treat 
the 0{p^) leg's as free parameters). 
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From the above simple considerations follows that in order to avoid potential problems with 
the instabilities of the chiral expansion, which might be present in the three-flavor xPT^ i^i the case 
of small X(3) and Z(3) (cf. (|lj, [ij)), we should g | 



• carefully define the bare expansion 

• confine ourselves (as far as the bare chiral expansion is concerned) to the linear space of "good" 
observables and keep the "dangerous" observables in the nonperturbative "resummed" form 

• use rather S(3), fb(3) (or X(3) and Z(3)) and r = 2ms/ {mu + md) as free parameters^ 
instead of expressing them in the form of the series in PGB masses and decay constants 

• eliminate the 0(p'*) LEC's algebraically, using bare expansions of "good" observables such as 

In the next section we shall illustrate the possible subtleties of the first step of this general 
recipe on the concrete example® of the PGB scattering amplitude P1P2 ^ PzP^ 

3. Construction of the bare expansion for the scattering P1P2 — > P3P4 
3.1 Chiral expansion of the "good" observable 

Let us assume a scattering of pseudoscalar mesons P1P2 P3P4 with masses Mp.. The amplitude 
S'(s, t; u) is defined as 

{P3{k3)P4iki)out\Piiki)P2ik2U) = i(27r)4<5(4)(A:3 + ^ h ~ k2)S{.s,t;u), (3.1) 

where s, t and u are the usual Mandelstam variables. The amplitude is related to the "good" 
observable^ 

4 

G{s,t-u)^\{Fp^S{s,t-u), (3.2) 

1=1 

(where Fp. are the decay constants) which can be directly obtained from the (cut) four-point 
function of the axial currents. Let us write for G{s,t;u) the following strict chiral expansion in 
terms of the low energy constants 

G = G(^) + g'^ + g|^', + Git\, + GSg. (3.3) 

GSg accommodates the higher order remainders. Using the functional method, G can be obtained 
from the generating functional 



F^Z[U,v,p,a,s] = F^ J d^x (c^'^\U,v,p,a,s) + C^^\U,v,p,a,s)^ 

+F^z[2p[U,v,p,a,s] + ... (3.4) 



^Note, that r is related to the "dangerous" observable 

2 - l = r + ... 

^We will do it for L4 — Lg but leave Li — L3 free, also L7 is a special case, see in what follows. 
*We shall tacitly assume the case of three light flavours in what follows. 

^Strictly speaking, the "good" observables correspond to the subthreshold parameters derived form G{t,s;u) in 
an unphysical point away from singularities. 
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by setting v — s=p = 0, s = 2BoA4 and expanding in the fields $ where U = exp(i$/i^o)- 
Following the notation in j|] , we have 

Ziolp[U,v,p,a,s] = z'^2i[U,v,p,a,s] + zl^l^^[U,v,p,a,s] 

= Indet + ~ '-Tt{D^^5D^'5) + . . . . (3.5) 

In the above formulae, 

D'^^ = 6''^n + ^Bfitr{{\'',\^}M) (3.6) 

and Ai is the quark mass matrix. Note this representation of Z^'^'l^ assumes that the masses running 
in the loops are the 0{p^) masses rather than the physical masses. Or, in more detail, provided we 
start with the chiral expansion of the squared product of the masses and decay constants 

i^|A/|, = (F|Af|,)(^) + {F^Mlf'^^ + F^MISfm^, (3.7) 

the masses in the loops are defined as 



2 (F2M|M2) 



Mp^ ' \r ■ (3.8) 





Note, however, that this is the first term in a potentially "dangerous" expansion of the ratio 



FlMl [FlMlY^) + (f 2M2)(4) + FlMlSp _^ 



From this definition of Z^^^,^ we obtain G^''-' = G^j' + G|^|; + G^'^^it which is exactly renormalization 
scale independent even for the external momenta off-shell. This meets the requirement of the 
renormalization scale independence of the remainder 5g- 

The first two terms of the above strict chiral expansion for G{s,t;u) have a serious drawback 
in the sense that the singularities in the complex stu planes required by unitarity are not placed 
at the physical thresholds but rather at points given by the leading order terms Mp of the chiral 

o 

expansion of the PGB masses. Straightforward substitution Mp^ Mp in the propagators of the 
loops, which apparently means merely a redefinition of the remainder Jg, could, however, in general 
spoil its exact renormalization scale independence. It is therefore desirable to use the freedom in 
the definition of the remainder more carefully in order to reconcile both scale independence of G*-^-* 
and unitarity. For this purpose, a useful tool is the matching with a dispersive representation ||^ 
of the amplitude 5(5, t; u) based on the reconstruction theorem |29[| . 

3.2 Dispersive representation for G{s. t; u) 

The above mentioned reconstruction theorem for the PGB scattering amplitude is based on the 
basic properties of unitarity, analyticity and crossing symmetry and provides us with the most 
general form of the PGB scattering amplitude up to the order 0{p^) in terms of dispersive integrals 
with known discontinuities. It was first proved for the case of tttt scattering in and for ttK 

scattering in |3^, ^ and since then it has been intensively used in various contexts. Here we use 
the general form of the theorem, more detailed discussion of which will be presented elsewhere . 

For the scattering of pseudoscalar mesons P1P2 — * P3P4, let us denote the s— , t— and w— channel 
amplitudes as S{s,t]u), T(s,t;u) and U{s,t;u) and write their partial wave expansion as 

00 

A{s,t;u) = 327r^(2/ + l)ylz(s)P; (cos 6*^), (3.10) 
1=0 
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where S, T, U and 

^=0^^^ = — T7J7T7T7J7— ■ (3.11) 

Here Ai(s) are the partial waves, 

Xa^Js) = {s- {Mp^ + MpJ'){s - {Mp^ - MpJ') (3.12) 

is the triangle function which corresponds to the initial/final state Aij (consisting of the pseu- 
doscalars PjPk) of the process in the channel A and 

^Mj=Ml^-Ml^- (3.13) 

According to the theorem, we get the following representation for the amplitude S{s,t;u) 

Sis, t- u) = S{s, t- u) + S^mt{s, t- u) + 0(/), (3.14) 

where S{s,t;u) is a third order polynomial with the same symmetries as the whole amplitude 
S{s,t;u). The nontrivial analytical properties are incorporated in the unitarity part iS„„it(s, u), 
which can be expressed as 

Sumt{s,t;u) = $^(s) + $^(i) + $^(u) 

+ [s(i-u) + Ai2A34]^'^(s) 

+ [t(s-M) + Ai3A24]*^(i) 

+ [w(<-s)+Ai4A23]«'^(m). (3.15) 

In the last expression, Aij = - Af|,.. The functions $-^(3) and *^(s) with A = S,T,U are 
analytic in the cut complex plane with the right hand cut from ta = mini j(Mp. + ^^PjY (where 
PiPj are the possible intermediate states in the given channel A) to infinity with discontinuities 
given by the formulae 

disc $-^(5) = 327r6l(s-rA)discylo(s) (3.16) 
disc^''^(s) = m^eis - TA)Aisc--pr^^^^^ — . (3.17) 

Ayf(s)Ay;(5) 

Here Ao(s), ^i(s) are the corresponding 1 — 0,1 partial waves. 

Consequently, once the right hand sides of (3.16, |3TT7| ) are known, the unitarity part Sunit{s, t; u) 
of the amplitude can be uniquely reconstructed to 0{p^) up to the polynomial, which encompass 
subtraction polynomials for the dispersion integrals. 

Let us now assume the chiral expansion of the amplitudes in the form 

A{s,t-u) ^ A'~'^\s,t;u)+ A''^\s,t;u) + A5a, (3.18) 

00 

A(")(s,t;u) = 327r ^(2/ + 1)A|"^(s)F, (cos 61, (i)). (3.19) 

(=0 

Starting from the 0{p^) amplitudes, we can use the two particle partial wave unitarity to get the 
discontinuity of the partial waves A'f'' (s) along the right hand cut^*' 

discA[^)(s) = ^ A^Jf(flAp)«^-4/(^)^(2).,^A.(^)* ^ ^^^6^ (3 



^It can be shown that more than two particle intermediate states yield contribution of the order 0(p*) and higher. 
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Here Zij = 1, 2 is a symmetry factor taking into account the possibility of identical particles in the 
intermediate state ij. Inserting this into the dispersive integrals we easily^ ^ get a minimal form for 
the 0{p^) unitarity corrections in terms of the functions $('*''^(s) and "^^^^^{s) reconstructed from 
the 0{p^) amplitudes 



Provided A, 



X|/(4)A(, 



or 



^327TfY,—J,,{s)A, 



1 ^ . . .i2)A.^^,^^^^(2)^,^A,^^y 



(967r)2 



1 = , 



>^A-{s)\'i^is) 



(3.21) 
(3.22) 



0, which will be our case, it can be shown that we only need one 
subtraction, Jij{s) = JlAs) — JLiO) instead of Jij{s). The explicit form of the function Jy (s) is 



given in the Appendix 

The above formulae can be used to write a dispersive representation of the "good" observable 
G(s, t; u) — Jli=i PPiS{s, t; u) to the next-to-leading order in the form 

G{s,t;u) = g{s,t;u) + Gunitis^t^u), (3.23) 

where Q{s,t]u) is the polynomial part and the unitarity corrections up to 0(j)^) are included in 



+ [.s(t-u) + Ai2A23]V'^(s) 
+ [t{s -u) + Ai3A24]^^(t) 
+ [u{t - s) + Ai4A23]lp^ (u). 



(3.24) 



Our goal is to write down a representation of t/i^^)^ and ip'^'^^^, which, notice, are distinct quantities 
from and ^S^^^^, analo gous to ( [3.2l| , |3.22D . Note, however, that while the relation of G(s, t; u) 

and S{s,t;u) is unambiguously fixed to all orders by (3.2), the amplitude can be defined order by 
order in various ways. For example, for the "good" observable G, the leading order piece G*^^^ of its 
strict chiral expansion is fixed by the lowest order Lagrangian C^'^\ but the corresponding 0{p^) 
piece of the amplitude S can be related in various ways. Similarly, the same is true order by order, 
where the amplitude at the given order can be defined up to higher order corrections. 

The most straightforward way is to write a safe expansion for 5(s, t; u) in the form 



S(s, t; u) 



^4 \ ~1 

Y[ Fp^ (g(2) (s, t- u) + G(4) (s, t- u) + GSc 

\i=l ) 



(3.25) 



with physical values of Fp. , thus satisfying the relation (3.2) order by order 

5(")(s,i;w)= (n^^^l G(")(s,i;w). 



(3.26) 



As we'll see, the minimal modification of the form derived from the generating functional is obtained 



l^Note that are real polynomials of the first order in s, t and u. 



- 11 - 



by using an alternative, potentially "dangerous" expansion 



S{s, t; u) 



G{s, t; u) 



(n ^o(l + + ...)) (g(^) is, t, u) + t;u) + ..) 

F,^G^'^ is, t; u) - iF„-*^G(2) t; u) ^ [F^p^ ) ^ + F,-^G^^^ [s, t; u) + 



(3.27) 



which defines 



S^'^\s,t-u) = F^'^G'^'^\s,t-u) 



(3.28) 
(3.29) 



i=l 



The representation of (f)'^'^'^^ and i/'^'*'"^ is therefore not unique. According to our definitions of 
the amplitude we get either 



1 Jij{s) {2)Ai~,i3 1 .^(2)ij^Ai 



p2 F2-G^ 



{s)G\; 



is)* 



■^11 J- p.-*- p. 

\2 1 T../'o^ ^(2)Ai^ij, 



(4)^ _ (96^)^^ 1 Ms) Gr-'-^(.)Gp-^^(.)* 

^ 2. ...^2^, Ay/(.)Ay;(.) ^ 






(3.30) 
(3.31) 



corresponding to the definition (|3.26|) or 



1 



^ (327r)2F-4^^J,,(.)G(^'-^--^'^(s)G|,')'^--^^(s)* 



(967r)2^_4Y- It ,.,GP'^--'^(.)Gr---^(.)* 



,1/2. N,l/2, , 



(3.32) 
(3.33) 



when reconstructing the bare expansion of Gfrom the "dangerous" expansion (3.27) and using the 
definitions ( |3.2S| , |3.29| ) for the 0{p'^) and O(p^) amplitudes. 

3.3 Matching the strict chiral expansion to the dispersive representation 



The dispersive representation (3.23) can be now matched to the formula (3.3). As we have mentioned 
above, the positions of the cuts in the formulas (3_^) and ( 3.23| ) are not the same; in the former case 
they correspond to the 0{p^) masses (3.S), which ensures the renornialization scale independence, 
while in the latter they are determined by the physical ones, as required by the unitarity conditions. 
In order to reconcile both these requirements, one can proceed as follows (c.f. also [||). 

In (^^), the nonanalytic terms are generally of the form P{s)J^j{s), where Jlj{s) is the renor- 
malized scalar bubble defined in Appendix ^ and P{s) is some second order polynomial. As the 

first step, one rewrites these expressions in terms of Jij (s) writing Jlj (s) = JJ^ (0) + s (0) + Jij{s). 
This adjustment allows us to split G uniquely into a polynomial part Gpoi and a nonanalytic part 
Gcut which accumulates the unitarity cuts 



G{s,t;u) = Gpoi{s,t;u) + Gcut{s,t]u) + GSq, 



(3.34) 
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where 



Gpoiis,t;u) = (G(s,t;u) - GSg) Ij^^jj^^^^- (3-35) 



Both parts are now renormahzation scale independent. 

As a second step, we replace the Gcut{s,t]u) with Gunit{s,t]u) from ( 3.25 ). This means we 
write 

G(s, t] u) = Gpoi (s, <; u) + Qunit (s, i; m) + G5q (3.36) 

where (5^ is a new remainder defined by this equation. According to the naive chiral power counting, 

Gcut{s,t;u) - Gunit{s,t;u) ^0{p^). 

The third step, not necessary from the point of view of preserving unitarity and renormahzation 

scale invariance, consists of a further modification of Gpoi{s,t]u) by means of replacement of the 
o 2 

0{p^) masses Alp in Jlj{0) with the physical masses Mp. This replacement does not spoil the 
renormahzation scale independence of the Gpoi{s, t; u) and corresponds to the convention introduced 
in This again means a redefinition of the remainders Sq, i.e. re-shuffling of the terms of the 

next-to-next-to-leading order. 



Note that the origin of Jj'^ (0)'s in one loop generating functional (3^) is twofold: they can stem 
either from the tadpole part Z^^^^ or from the unitarity corrections Z^j^^u ■ argued in |^ that in 

the former case the above mentioned replacement does not necessarily modify the numerical value 
of the remainders much. The reason should be that the chiral logs appear only in the combination 

o 2 o 2 

lip (xAIp ln(Afp )• The replacement here means 

a/p ln(Mp /m^) ->A/p ln(Af|,//x2). (3.37) 

o 2 o ^ 

Because Alpoc Y, the difference should therefore either be small for F ^ 1 (where AIp^ Mj,) or 

the contribution of /ip itself is tiny for Y ^ 0. 

On the other hand, the logs from Z^^^^^ do not generally come with such a prefactor. Therefore, 
o 2 

with a replacement A/p— + Mp inside Jlj{0), one might create large differences between the "old" 
and "new" remainders due to the enhancement of the contributions of chiral logs for small Y — X/Z. 
However, without the replacement inside the chiral logs of this type we could expect an unphysical 
increase (and irregularities) of the observables for Y ^ 0. Also, here the replacement is natural 
physically, remember that the matching with the dispersive representation consist essentially of 
an analogous replacement within the unitary corrections. Let us also note that the splitting of 
the generating functional into the tadpole and unitarity part is not unique (it depends e.g. on 
the parametrization of the fluctuations around the classical solution of the 0{p'^) field equations 
in the functional integral), though the sum must be independent on this and therefore it is more 
consistent to use the same rule for both^^. Nevertheless, it could be of some worth to test the 



differences between various treatments of the chiral logs numerically (see Subsection S.2 ) 



The resulting bare expansion ( 3.36 ) now not only meets the requirement of the exact scale 
independence of the remainder Sq, it has also correct physical location of the unitarity cuts. Of 
course, we could achieve the last property simply by inserting physical masses into the functions 
Jij(s) in Gcut{s, t] u). The replacement Gcut{s, t; u) — > Gunitis, t; u) has however another advantage. 



Namely, using the prescription (3.3C), ( 3.31| ), the corresponding amplitude, written in the form 



^■^Also notice that the offending Y dependence of the chiral logs with O(p^) masses inside comes always in the 
combination Y I ^ where is the renormahzation scale. Provided we were able to reparametrize the bare expansion 
in such a way that all the running 0(p*) constants were completely expressed in terms of the physical observables, 
the explicit independence on [i would at the same time guarantee an elimination of the irregularities for y — » 0. Such 
a treatment has to include the reparametrization of L1-L3, which is, however, beyond the scope of our paper. 
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(without any expansion of the denominator) 



S[s,t;u)^—^ — — (3.38) 

lU=i 



satisfies the relations of perturbative unitarity (with S'^'^^ and S''^^-' given by (3.26)) 



disc5,(^)(.) ^ —^^^s[^^''-^^{s)s[''^'^-''\sr (3.39) 

ZiT s 



exactly (i.e. not only modulo the next-to-next-to- leading correction), which can be sometimes 
technically useful {e.g. for the unitarization by means of the inverse amplitude method, cf. pi). 
The same is true using the prescription (B.32), ( ^.3S ) with S'^'^^ and S^^^ given by (3.28, ^.29] ). 



As we shall see in what follows, the latter prescription gives a minimal modification of the strict 
expansion ( |3.3| ) compatible with exact perturbative unitarity. 

4. General properties of Tzrj scattering amplitude 
4.1 Basic notation 

Let us denote the s— and u— channel amplitude in the isospin conservation limit as 

{^'{p')v{q)outU^{p^)fl{p\n) = i(27r)4^(F/ - F,)<5"^5(s,t;u) (4.1) 

and the crossed amplitude in the channel as 

{ri{p)r){q)out\n%p'')n\p%r.) = i(2^)^5(F/ - Pi)5^'T{s,t-u). (4.2) 

Crossing and Bose symmetries then yield 

T(s, t; u) — S{t, s; u) 
S{s, t] u) — S{u, t; s) 

T{s,t;u) = T{s,u;t). (4.3) 

Writing the partial wave expansion as 



S{s,t;u) = 3271^(21 + l)Pi{coses)Si{s 



1=0 



(t - u)s + 

the scattering lengths a/ and phase shifts 6i{s) are given by the formulae 

ReSi{s) = ^P'^\ai + (^(P^)) for P ^ , s ^ (M^ + M^f 
fAP \ 

5i{s) = arctan ^ReS';(s) , (4.5) 



where P = A^(-^(s)/2y^ is the CMS momentum. I.e., in the units of (pion Compton wavelength)^'^-'^ 



ai = Af2'+i hm ^ ReAi{s). (4.6) 
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Let us also define the subthreshold parameters Cy in terms of the expansion of the amplitude in 
the point of analyticity t = 0, s = u = E^,r = + 

= ^c,jtV2i (4 7) 

where 

"=4M7 = AM, ■ (^-^^ 

The dimension Cy is dim[cy] = raass~'^^~'^'^ , in what follows we will refer to the dimensionless 
numbers CijM^^'^'^^ . Let us note that in the limit rriu = rud = we have two Adler zeros at p° = 
and = 0, which implies the following SU{2)l x SU{2)b. theorem 

lim Coo = 0. (4.9) 
We can also quote the low-energy current algebra result 

Sis,t;u) = ^, (4.10) 



which is in agreement with (4.9) 



4.2 Dispersive representation 

As a result of the symmetry properties of the amplitudes, the dispersive representation to the 
next-to-leading order ( |3.24 ) for GT,,(s,i;M) — F^ F^S{ s, t; u) simplifies, namely cj)^ = (p^ = (j) and 



= = ip. The intermediate states in ( |3.30 , 3.31 ) are^^ ttj] and KK in the s and u channels 



and TTTT, rji] and KK in the channel. This implies V'(s) = 0{p^), because the P— waves in the 
s— channel start at O(p^) due to the low-energy theorem ( 4.10 ) for the irrj — > irrj amplitude and as a 



result of charge conjugation invariance of the ttt] —> KK amplitude. Moreover, tp'^ = 0, because the 
partial wave decomposition of the channel amplitude T{s,t;u) contains only even partial waves 
due to Bose symmetry and charge conjugation. We therefore get 

6*^^(5, t;u) = GT,n..poi[s,t\u) + G7vri,unit{s,t;u) + 0{p^), (4-11) 

where the polynomial part has the following general form 

G^^,poi{s,t;u) =a + iSt + '^t'^ +uj{s-uf. (4.12) 

Note that the parameters a, . . . ,lu are related to the expansion of the Green function G-j^riis, t; u) 
at the point of analyticity t — 0, s = u = S^^r = M!^ + and therefore they represent "good 
observables" according to our classification. 
The dispersive part is 

G^r,.umt{s, t- U) = (t)^{t) + (f>{s) + <j){u), (4.13) 



where (t)^{t) and 0(s) are given by the formula (3.30). A complete list of relevant leading order 
contributions Gq^^^^^'-' and Gq^)*"'^^'' can be found in Appendix |c[ here we give the resulting 



'Here we assume isospin conservation. 
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expressions (transcription to the convention (p.32|, 3.33) is straightforward) 



glV 3 3 3 

-(2 Mk - - Mr,)? ^ ' 



1 o2 / o2 o2\ J (s) 



X [(3. - 2M| - 2Af,^) + (2 A/' mL)]^^^^} • (4.14) 



K 



In terms of these functions, we have (notice that (/> (0) = 0) 

^ 2j.2 (M.^V,) + 16^^'^' + '^((^^- + ^'^)') + -^((^-^ - M^f)) 
1 Af'^ / / ' \ 



and 



^=00 = Y2WI + 20(S^^)) 

?7 TT 



77 TT 

= ;^(7 + ^'^'^"(0) + ^/(S„.) 
1 6M^ f 1 

r2 p2 
?7 7r 



COI = ;^(^+4</> (S..)). (4.16) 



While the scattering lengths, being related to the value of the amplitude at the threshold, are 
not candidates for "good observables" , the situation is a little bit more subtle in the case of the 
subthreshold parameters. Provided the rj decay constant was known from experiments as accurately 
as i^TT, then (similarly to a, . . .) also the Cy could be treated as "good observables". However, this 
is not the case, and we should rather use a chiral expansion of F^i in the above formulae. Therefore, 
the subthreshold parameters are typical examples of the dangerous ratios, which should be treated 
with care. 

4.3 Bare expansion for G{s,t;u) 

For the strict expansion in terms of LEC's (i.e. without any reparametrization in terms of physical 



observables) derived from ( 3^ , ^^) , we have confirmed the results of the article ||3^ by independent 
calculation. The 0{p'^) expansion can be written in the form 

G., - G(2) + Gi? + gI^), + GitL + GSg, (4.17) 
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where 



G('\s,t;u) ^ ^ All 



3 

G'-^\s,t;u) = 8{L[{^,) + ^Ll{^,))it-2M'^){t~2M^) 

+4(L^(m) + lLli^^ms - Ml - M^)^ + _ Ml - Ml^f] 
+m^{Mt - 2Ml) M, +{t - 2M2) mI] 

g Q 2 Q 2 Q 2 Q 2 

--Ll{fi){Ml + Ml,) M^ +8i^(M) AU (M. +5 M J 

o2 o2o2g4 o4 

(41 1 o 4 

+ - Ml - Af2 + ^ m!]V],^(s) + ^[zi - Af2 _ Af2 + 2 ^^J2 
+ 1^! [t-2Af2 + ^M!]j;jt) 
+ - (A/^ Mjj;,(i) 

+ - 2Af2 + 2 A/J[3< - 6Af2 + 4 A/'] J],kW (4.18) 

are the 0(p^), counterterm, tadpole and unitarity contributions respectively. In the above formulae, 
the masses within the loop functions Jpqit) are the 0{p^) masses 

o 2 o 2 o 2 2 

M^=2Bom, MK^B^m{r + l), M^= -Bom {2r + 1) . (4.19) 
The chiral logs fip can be expressed using Jpp(O) 

Written in such a form, the sum G^^'' +G'|^);+G'|j'[]jf is exactly renormalization scheme independent by 
construction. Let us now proceed as described in the previous section and write the bare expansion 
of G(s, t; u) in the form 

ttt;, pol 

{s,t;u) + g {s,t;u) +G6a. (4.21) 



Writing J[j{s) — Jij{0) + Jijis) in (4.18), we get the renormalization scale independent polynomial 
part 

G^^,poz(s,t;u) = G^-^\s,t;u) + G'^^\s,t;u) + GfJ^is,t;u) 
+ ^mI [t~2All + ^Ml]JU0) 



2 o i 2o2o2 1o2 

Q AT, J;^(0) + - (A/^ -- AU] 
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-[t - 2Ml + 2 Aif!][t - 2Ml + ^ m!]}J;,k(0). (4.22) 



Comparing this with the general form ( 4.12 ) of Gnri.poi{s, t; u), we get for the bare expansions 
of the parameters a — lo the following manifestly renormalization scale independent form 

l„o2 1 o2/7o2 11 o2 



" = 3^° "^'^ +96^ U ^^^^ 

+4[8(LUa*) + ^i5(A*)) + ^ J;,k(0)]M,2A£,^ 

-[16i^(M) + J],K(0)]Af^ Af' -[16i^(M) + Imf,) + ^ J],k(0)]M2 
-llmfi) - i J],^(0) + ^ J;.(0)]M2 mI 

+[40ig(M) + ^j;,(o) + ^j^KKm mImI 

2 2 2 

+32L^(m) {AU - M^) 

+[^Li{n) + yL^(M) + j;.(o) + ^j;,(o) - ^j;,(o) + ^j],^(o)] mI +IfImX 

(4.23) 

+ [8LI{p){mI + mI) + IrKKiO) Mr, +Uj:AO) + ^JkAO)) Ml]+Ph (4.24) 



7 = [8(LUm) + \Ll{^l)) + \jWm + [2(i^(M) + Ji^(A*)) + ^JWrn + 75, (4.25) 
[2{Ll{^i) + iL5(^)) + ^ J;,^(0)] + (4.26) 



5. Reparametrization of the bare expansion 

Let us now discuss the various possibilities of the reparametrization of the bare expansion. 

5.1 TTTy scattering within the standard chiral perturbation theory to Oip^) 

The standard way of dealing with the chiral expansion consists of two "dangerous" steps. The first 
one involves using the inverted expansions of the type ( 2.10 ) in order to express the amplitude in 
terms of the masses and decay constants instead of the parameters Bofh, Fq and r — ms/fh of the 
Oip^) chiral Lagrangian. Here one encounters an ambiguity connected with different possibilities 
how to choose the observable G in (2.10), the chiral expansion of which starts with the desired 
O(p^) parameter Gq. 

Let us fix this ambiguity by using the expansions of F^, M^ and M^, inverting of which leads 

toi4 

= (1 + 4^^ + 2/iK) - (LK/i)(2 + r) + Ll{p)) (5.1) 



* Instead of we could use the chiral expansion of A/^ to obtain 

_ 3 / M 2 1 
'■ = '■^ = 1 Mf-3 
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2Bom = Ml 



r2 



+ \^^'n ^(2i8(M) + 2(2 + r)L^(M) - ^^(m) - (2 + r)UM] 



9 Af2 



(5.2) 
(5.3) 



Inserting the inverted expansions (5.1-5.3) into ( 4.23 ) and ( 4.24 ) and keeping terms up to the order 
0{p'^) we get 



1 



1ft «4 
-Fl;Ml + -MlMl^Llii,) - -U,{i,)Mt{r2 - 1) 

+MlMl{Z2L\[i,) - 16L^(m) - \Ll{ii)\ 

+ i(2r2 + \)M%2Ll{^) - 16iK/i) + ^j;^(0)] 

+Af;f[-^L^(^) + i6l^(m) - \JIM + ^^;,(o) - ^ j;,(o) + \j],Km + c'^; 



(5.4) 



1 



M2[16L^(A.)(r2 + 2) + J]^^(0)(r2 + 1) + JIM\ + /^^ 



'/3 ' 



(5.5) 



with new remainders (5** and (5^* 



which might be, however, out of control as we have already 
discussed. In fact, this first step involves three "unsafe" manipulations from the point of view of 
resummed xPT: using "dangerous" expansions for the masses as a starting point, the inversion 
and finally the negligence of all higher order terms generated by this procedure after the insertion. 

It's understood to use physical masses inside the chiral logarithms. Higher order LEC's are 
then fitted by using additional experimental input, no parameters are therefore left free. Also note 
that (^.3[) effectively implements the classical Gell-Mann-Okubo formula 



3M2 - 4Af 2 



K 



Mi - 0. 



(5.6) 

This insures renormalization scale independence. We, however, leave at its physical value in 
cases when it was produced by on-shell mass on outer legs or inside chiral logarithms, which is 
compatible with the requirement of scale independence. 

The second step is connected to the fact that the amplitude is used in standard xPT^ rather 
than G{s,t;u). As was shown in Section 3.2, the expansion of the amplitude can be organized in 
various ways, of which only ( 3.26D is considered safe in the resummed approach. On the other hand, 
from the standard point of view it often seems more advantageous to use ( 3.2^ , 3.2S| ), as together 
with (5.1) it leads to only the experimentally very well known pion decay constant being present 



in the formulae. This can be seen on the case of -F"^, which is experimentally poorly known due 



to rj — rj' mixing |43 and thus if it's kept at its physical value as was done in ||38| , a significant 



uncertainty is introduced into the results. As the normalization ( 3.2S , 3.29 ) is used more often in 
NLO SxPT, we will adhere to this view and perform this second step by expanding the kaon and 
eta decay constants from the denominators and subsequently cutting off the higher orders. 



or even F^Ai'^ to got 



The latter choice, formally as good as the previous two, could also involve the redefinition of the loop masses to 

Mp= FpMp/ instead of the simple Mp= Mp as in the case of the other standard reparametrizations. Even then, 
however, it suffers from numerically large 0(p'*) corrections which could produce instabilities of the reparametrization 
based on this observable. 
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Using therefore the prescription ( 3.32 , 3.33 ), the dispersive part of the 0{p'^) amphtude ( 4.14 ) 

ti 



simphfies using the reparametrization recipe described above 

1 



9 

+ '-M'^ - Ia/J - ^A4) - lMUr2 iriKKis), 

0^(s) = ^M^{s - ^M^)X4s) + ^M* (8r2 + l)J,,(s) 
+ - l^Ml - ^a4) + \Ml(r2 + 3)] 

X [(3s - 2Af|. - 2Af2) + iA/2(3r2 + 1)]Jkk{s). (5.7) 

The second step also propagates itself to the case of the subthreshold parameters and the 
scattering lengths aj, where it consists of the expansion of in the denominator of the formulae 
( 4.16D and (4.15). This step could in principle produce uncontrollable contribution to the remainders 



as well. 

5.2 Resummation of the vacuum fluctuation 

In order to preserve the global convergence, as was discussed, in the context of resummed xPT we 
are not allowed to perform "dangerous" inverted expansions and thus to express the 0{p^) masses 

o 

M p and the decay constant Fq in terms of the physical ones in the way it is common within the 
standard xPT calculations sketched above. Instead of this, the O(p^) LEC's are left free, or more 
precisely, rewritten using parameters directly related to the order parameters of the chiral symmetry 
breaking^^ 

m, _ F§ Fg 

The bare expansions for masses FpMp and decay constants Fp are used the reparamctrize the NLO 
LEC's L4 — Lg. As the dependence is linear, it can be done in a purely non-perturbative algebraic 
way by introduction of an unknown higher order remainder to each observable used. The relevant 
formulae for L4 — Lg can be found in Appendix ^. 

As masses and decay constants do not depend on Li, L2, L3, bare expansions of some additional, 
experimentally well known observables is needed for these LEC's. This is, however, even if highly 
desirable, out of the scope of our article and we make a shortcut and use the standard tabular 
values for these constants. We will make an analysis of the sensitivity of our results to a change in 
the value of Li — L3 in the next section devoted to numerical results. 

For the resulting expression for the parameters a and /3, we use the following abbreviation for 
some repeatedly occurring combinations 

P2 ]\,f2 

sir) = 2%^ (5.10) 
— 1 



= ^ ( ^1 - 1 ) (5.11) 



p2 

^ TT 

Agmo - ' ' ;2M2 (5-12) 



^Here we omit the explicit dependence of X, Y and Z on Nf keeping in mind that Nf = 3 in what follows. 
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in terms of which we get 

+2^:^(2 - 1) (jl'il^'- + 1) + ^'n) + jf2Vir) (vMl - i(r - 4)M; 

1 X -. , , , 3 X , ,9 , --9 1 f X 



967r2 Z 

Mi X 



M^(Ar + 5) + -^—M^M'^ — ( M'^(4Ar + 67) 

[j;^(0)(2r + 1) + 2rj,,,{0)r - j;,(0)(4r + l)]r 
[j;,(0)(2r + l)(r ~ 4) + j;,(0)(19r - 4) - 2J],^(0)(r2 + 6r - 4)] 



2(r 


+ 2)(r- 


1)Z 




M2m2 


X 


6(r 


+ 2)(r- 


l)Z 






-( 


18(r + 2)(r 





^X2 



- j [J;,(0)(5r2 - lOr - 4) + 6J],^(0)(3r2 - 2r - 4) 
+4j;„(0)(r2 + r - 2) - 9j;,(0)(3r2 - 2r - 4)] 

(5.13) 

P = |F2(l-Z-77(r)) 

1 )i#2 V" 1 

+ 3 ^ YK,(0)(2r + 1) + JWmr + 1) - J;.(0)(3r + 2) - ^(r + 2)(r - 1)] 

-2E,.[8(LJ(Ai) + \Ll{^i)) + ^ J;,^(0)] + pS'p. (5.14) 
The new primed remainders are the following functions of the original remainders entering the game 

2 / Af2\ p2 

^3r-(r-4)-| (^^,,-<5,J 



(r + 2)(r- 1) \^ ' 'M^J 'F^ 
2r(2r + 1) / ^2 M2 ^ , 

- (. + 2)(.2_1) (2^^..-M.-(^ + l)^^.M. 



3t^'5^.m. + Sf^m^ - ^^I^^F^M^ 1 (5.15) 
a' a , 2 2f 2.J^, - (r + l)F2^j^^ 

This is an alternative to the first step used in the standard approach to xPT . Because the value of i^^ 
is not very well known, we make in a sense a parallel to the second step as well, i.e. using the chiral 
expansion for F'^ in the denominator of ( 4.14| ), ( 4.15 ) and ( 4.16 ). It involves the reparametrization 
in terms of X, Z, r (see Appendix]^ for details), but, contrary to the standard case, the denominator 
is not further expanded and the result is given in a nonperturbative resummed form of a ratio of 
two "safe" expansions. 

5.3 TIT] scattering vifithin the generalized chiral perturbation theory to 0(p*) - the bare 
expansion of G(s, t\ u) 

In analogy with ( [4.18|) , the strict chiral expansion for G{s^ t; u) within the generalized can be 
straightforwardly obtained by using the Lagrangian summarized in the Appendix ^, where we use 
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the traditional notation for the LEC's. The result has the following structure 



where^^ 



G^^\s,t-u) = -F^[Ml+Am^{:iA,-^r-l)ZP + 2{2r + l)Zl)] 
G^^\s,t;u) = ^F^ [-2m (BM^ + M^ (2 + 4r)-2 (2 + r) t) |-2mE^^e 



+ 81 rh pi + m p2 + 



64 r- 16 r^) 



+ (100 + 64 r + 34 r^) fh^ + (2 + r^) 

+ (96 - 96 r) pe + (144 + 288 r + 108 r^) fh^ pr] 



(5.17) 



G';^\s,t;u) 



8{Li + -L3){t - 2M^){t - 2M^) 



+4(i2 + -L3)[(s 

+\fn?Fi {-(Si - B2)T.^r> + 2DPMl{r - 1) - 2CfM^{r - 1) 

+ Cf (2r + l)t - D^[]^^„i{br + 4) - (2r + l)t] 

- 2B4[3M2 + Ml{2r^ + 1) - (r^ + 2)t]} 

+ \fn^F^ [2b6Ei + 16£;2 + (256 - 256r2) + (32 + IGr^) 



J^f (256 + 320r2) + ^5^^(192 - 320r 

F2^(240- 216r-24r-'';TX6 
f^f (16 - 8r - 8r^) + ^1^(16 + lOr + lOr 
Fg^^(32 + 40r 



leOr^ - 32r3) 



+ ^6^^(32 - 32r + 16r2 - ler^) 



16H + 20r^) + i^f ^(384 - 160r - 256r^ + 32r^) 
+ F2^(400 + 234r + lAr^) + F5^^(576 + 720r + 480r2 + 168r^)] 

-^F^ [2mBo{nr, + &IJiK + 9Af7r) + 9>Aofh^{%iJir, + 3iJ,K{r + 8) + 48/i„) 
+ 4Z^m^ {p^{16 + 41r) + Mif(48 + 90r) + At^(96 + 45r)) 
- 16Z^m2(2/i^(5r - 2) + 3/iK(6r - 4) + 3M^(3r - 8)] 



[Ml + 4m2(3^o - 4(r - 1)Z^ + 2(2r + l)Zi)f[r^^{s) + j;,(«)] 



1 
9 

+- s 



2- 



-{r-l)fh\A, + 2Z^)Yfj,j,{s) 



+-[u -Ml- Ml, + -Ml - -{r - l)m\A^ + 2Z^)f rj,j,{u) 
+ '^[Ml + 4m2(3^o - 4(r - 1)Z^ + 2(2r + 1)^0^)] 
x[t- 2Ml + -^Ml + 10m2(Ao + 2Zl)]j:^{t) 
+^[Ml + 4m^3Ao - 4(r - 1)Z^ + 2{2r + 1)Z^)] 



i^AU the LEC's in the following formulae are the renormalized LEC's at scale /i. We have omitted explicit notation 
of this in order to simplify the expressions. 
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X - Ull + a2((8r2 + l)Ao + 8r(r - 1)Z„^ + 2(2r + Ij^Z^f)] 
+ i[t - 2M2 + 2Ml + 8(r + l)m2(Ao + 2Z^))] 
X [3t - 67\f2 + - 

+ ^(r + l)a2(3rAo + 2(r - \)Z^ + 2(2r + l)Zo^] J^^(f). (5.18) 

In the above formulae, the generalized 0{p^) masses (also present implicitly in the chiral logs fip 
and the loop functions Jpq{s)) are 

= 2[Bo + 2m(r + 2)Z^]fh + AA^ih? 
M| = [Bo + 2m{r + 2)Z^]fh (r + 1) + A^m^ [r + 1)^ 

^\[Bo + 2fh{r + 2)Zl]m{2r + 1) + ^Aafh^ {2r^ + l) + ^Zo^m2(r - if. (5.19) 
The unitarity part can be further split into the polynomial and dispersive part 

^unit ~ ^unit,pol ^unit.disp ~ ""tmit U''->./(0) ""nriit I J''^ J ■ (^u.z,u; 

According to the general prescription, the dispersive part can be replaced with that of the dispersive 
representation which has the general form 

5':lt=^^it) + ^is)+^iu), (5.21) 

where now 



2 ■ 

-lis '-M'^ \mI - \mI \[2Ml - Ml - M^^ + a^^KMl)?^-^ 



(5.22) 



+ - '^Ml - + '^{(Mk - M^f + 2a^KMKM^] 

X [3. - 2Mi - 2M'^ + a,K{2M'^ - p^^)]'!^^ 



(5.23) 



for (3.30, 3.31) and analogously for (3.32, 3.33). The coefficients a^r,, . . • parametrize the difference 
between the standard and the generalized cases, within the standard 0{p'^) chiral expansion their 
values are either one or zero. The dependence of these constants on the LEC's are given in the 
Appendix 

5.4 Observables of the ttt] scattering within GxPT to 0{p'^) - the reparametrization 

As it can be easily seen from the above formulae (in fact it is a consequence of the construction of 
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GxPT), after identifying the parameters of the Lagrangians, 



Bom 

^L;{^) - ^m^Ao (5.24) 

and defining the remainders using the physical masses inside the chiral logarithms and the loop 
functions Jpq, the generalized bare chiral expansions contain all the terms of the standard one. 
More precisely, the generalized O(p^) bare expansions include extra O(p^) terms, which are counted 
as 0(j)^) and 0(p^) within the standard chiral power counting scheme. As a consequence, after 
writing the generalized bare chiral expansion of a generic "good" observable g in the form 

g ^ g(2)MxPT _^ g[^),GxPT ^ g(4),GxPT ^ ggGxPT (5 25) 

and then collecting the "standard" terms together, this expansion can be formally rewritten as 

g = 5(2)^ + g(^)^SU ^ gS^^ gg^ ^ gg(G) ^ g^G^PT^ ^^^^^^ 



(G) ^ ^ 

mulated in gSg . In the case of the polynomial parameters a . . . lu (4.22-4.2(;), the two versions of 
the chiral expansion coincide for 7 and 

6^ = 5,«xPT ^ ^GxPT^ (5 27) 

while the "standard" remainders Sa, 5p can be split into an explicitly known part, which includes 
the extra "nonstandard" terms, and the unknown remainders inherent to GxPT 

= 5J°"^^{^,) + 3^^r(M) + (5.28) 

TV IT 

pSp = /3<5/3'°°^'^(/i) + m^Fi5'f^^{^l) + /3,5^'^^^. (5.29) 

Here the first terms correspond to the new loops and the second to the new counterterm contribu- 
tions. The explicit expressions for them can be easily extracted from the formulae of the previous 
subsection, the results are however rather lengthy and we postpone them to Appendix |^. 

Let us note that both S''°°p{^) and S^j^{pl) are generally renormalization scale dependent. How- 
ever, due to the running of the GxPT LEC's Aq, Zq , , ^ and ^ , which, after identification (5.24), 



where the identification (5.24) is assumed. The extra 0{p ) terms mentioned above are now accu 



is the same as in the standard case, the "standard" remainders Sa and S/j given by (5.28, 5.29) are 
/x-independent. Of course, the "true GxPT" remainders S^^^^ , . . . , S^^^'^ are scale independent 
by construction. That means the sum of the loop and counterterm contributions to the "standard" 
remainders is /i-independent too. 

The usual way how to handle the reparametrization of the GxPT bare expansions is quite 
similar to the standard one. The difference is that as there are three additional O(p^) LEC's, not 
all of them can be reparametrized using inverted mass and decay constant expansions. The solution 



'The reason is that they stem from the terms quadratic in the Mandelstam variables. 



- 24 - 



is to leave two of them free {e.g. r and ^ = /Aq). Consequently, the expansion is performed 
according to the generalized power counting scheme and the terms of order higher than 0(p*) are 
discarded. 

We will, however, not use this approach, but rather exploit the relation ( 5.26| ), i.e. sew the 



standard and generalized bare expansions together. The reparametrization is then an extension 
of the resummed one (Appendix ^), w here all the remainders of mass and decay constant bare 



expansions are split according to (5.26). We can use all the resummed formulae as they are exact 
algebraic identities, valid independently on the version of xPT . The generalized contributions to 
the remainders can be found in Appendices ^ and ||. The outcome for the parameters a and (3 is 
then obtained by simply inserting all the generalized results for the remainders (Appendices ^ 0, 
^ into the expression for 5^ and 5^ ( [S.15 , 5.16). 



After this procedure, the generalized LEC's are present only in the formulae for the standard 
remainders. Also note that 5a°°'P'^{^), 5 p^""'^'' [fjL) as well as the generalized loop contributions to 
mass and decay constant remainders depend explicitly on the 0{p^) LEC's Bq — XM^/2rh, Fq — 
ZF^, Aq, Zq and Zq .^^ So as the last step of the reparametrization, the remaining dependence of 
the generic "loop" remainders da°°P^{fj,), . . . ,etc. on the O(p^) LEC's Fq, Aq, Zq and Zq can be 
removed up to the order O(p^) using the leading order expressions 



fh'F^A^ ^ -^FlMle{r). (5.30) 

As a summery, our handling of the generalized bare expansion can be viewed in two ways - 
either as a partial estimate of the standard remainders present in the resummed approach or as 
a special treatment within the generalized framework, where the 0{p^) (and partly 0{p^)) LEC's 
are reparametrized algebraically at the leading order, while treated perturbatively at the 0{p'^) 
one. The numerical results including a simple estimate of the remaining NLO and NNLO LEC's 
are presented in Subsection 6.6, also Appendix ^ contains an illustrative example of applying this 
procedure on F^. 



6. Numerical results 



In this section we shall present the numerical analysis of the observables connected to the nrj— 
scattering amplitude and the results which illustrate the subtleties of the various versions of the 
chiral expansions described above. In the numerical estimates we use M^^ — 135MeV, = 
548MeV, Mk = 496MeV, ^ ^ Mp = 770MeV, F„ = 92.4MeV and Fk = llSMeV. For the 
calculation within the standard xPT, the 0{p'^) LEC's are taken from ref. [Q and |4^. In the 
alternative reparametrization schemes, where only Li, L2 and L3 remain among the free parameters 
and the other ©(p**) LEC's are expressed in terms of physical masses, decay constants and the 
indirect remainders, we again keep (though rather non-systematically) the values of Li, L2 and 
from the same references. The sensitivity on this LEC's might be then estimated by means of 
the variation around these central values. In this chapter we insert the physical masses into the 
functions Jpq(O) unless stated otherwise. 

^*More precisely, the loops depend on the "true O(p^) LEC's"ylo, Zq'^ {c{. Appendix the difference is however 
of the higher order in the generalized power counting. 
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6.1 The standard chiral perturbation theory 







103/3/M2 


10^7 


lO^w 
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1.68 


0.90 


-1.52 


2.24 
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45| 


1.91 


-0.68 


-0.23 


-5.03 
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2.48 


7.49 


3.31 


9.48 



This subsection discusses the predictions of 
the standard chiral expansion to the order 0(p*), 
which are summarized in Table ^ and |^. Let us 
start with the parameter a of the polynomial part 
of the amplitude. The relevant formulae from the 



Table 1: The standard 0(j3'') values of the poly- 
nomial parameters for the two sets of LEC's taken 
form references and In the last row, the 

sensitivity on the LEC's is (over) estimated adding 
the uncertainties associated with the LEC's jiij in 
quadrature (This is of course only a rough estimate, 
because in fact not all the uncertainties of the Li's 
are independent). 



Subsection 5.1 and the LEC's taken from 



19 



result in the following value 



0.683 + &t ) FlM-l 



(6.1) 



In this expression, the first term corresponds to 
the current algebra result a'-^^ = F'^M'^/i, while 
the second one represents the O^p"^) correction. 
The third term is the standard remainder, which might be out of control when X,Z and r far 
from r2, even if the bare expansion of a were globally convergent (let us remind that a is a "good" 
observable) as we have discussed in Section |[ Let us also notice the unusually large next-to-leading 
correction, which could also indirectly indicate the numerical importance of the remainder in this 
scheme. 

The actual numerical value of the NLO correction is very sensitive to a shift in the 0{p^) LEC's. 
The corresponding variation Aa is numerically 

= (3.38ALi -f O.56AL3 - 3.5OAL4 - O.3OAL5 -f 3.62AL6 - 3.42AL7 + O.lOALg) x 10^. (6.2) 



For example, using the 0{p^) analysis based LEC's from 
Table 1)20 

An- 

= 0.23. 



Aa 



instead of those from |Q, we get (of. 

(6.3) 



Note that the large coefficients in front of the L4 and Lq contributions indicate sensitivity of this 
observable to the vacuum fluctuations of ss pairs as mentioned in the Introduction. 





Coo /Coo 


lO^cio 


lO^cao 


lO^coi 




lO^ai 


1 
1 


44 
45 
A 


1 
1 


1.06 
1.51 
2.49 


0.91 
-0.67 
7.49 


-1.23 
0.07 
3.31 


8.27 
-3.36 
15.16 


1.96 
1.18 
3.21 


0.59 
-0.60 
2.80 



Table 2: The standard 0{p'^) values of the subthreshold and threshold parameters as in Tab. |l|. The 
parameters are given in their natural units described in the main text. 



Let us compare this case with the related "dangerous" observable, namely the subthreshold 



parameter cqq. From (4.16) we get 



Coo 



1 AP 

- (1 + 0.683 - 0.625 -I- 0.006) — f = 1.0640^0^^, 

3 i'-rr 



(6.4) 



where the individual terms are the leading order contribution — M^/3F^, the next-to- leading 



correction to the parameter a, the next-to-leading correction to induced by the expansion of 



^^This set of LEC's is used in numerical estimates unless stated otherwise. 

■^"Because the values of the LEC's Li based on the O(p^) fit include implicitly parts of the 0{p^') corrections, the 
large variation can be interpreted as a signal of the importance of the NNLO contributions to the parameter a. The 
same is true for other observables from the Table |l| 
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the denominator and the contribution stemming from the unitarity correction 0(s) respectively. 
The first two large corrections accidentally cancel here, this, however, does not automatically imply 
similar cancellation of the potentially large remainders (we have not written them down explicitly 
here). Also, the strong sensitivity of a to the variation of the LEC's propagates here giving 

$^ = ^-0-28ALaxl03 (6.5) 

and furthermore increases the uncertainty of the O(p^) correction. This strong sensitivity supports 
the possibility that the standard remainders for coo might be numerically larger than the next-to- 
leading correction. Namely using the LEC's from the 0(p®) fit which generates part of the 
O(p^) corrections to the reparametrized expansion of coo, we get 

^ = 0.45. (6.6) 

We can also check the sensitivity of the next-to-leading order contributions to the way we 
rewrite them in terms of the physical masses and decay constants (i.e. how we use the 0(j>^) 
relations generating here a difference of the order 0{p^)). Provided we insert the alternative 0{p^) 
expressions for r into the chiral expansions of a and cqo 

^.= 2(^^-1] =24.2 (6.7) 

r*2 = 2^|| - 1 = 39.4, (6.8) 
instead of the standard 0{p^) value r = 7-2 = 2M|./Af^ — 1 = 25.9, we get as a result 

5 = i(l + 0.601)F2M2 (6.9) 
1 

Coo = o(l + 0.031)-f (6.10) 



and 



3^ ' ■ ' 



a* = + 1.297) F^M^ (6.11) 
1 

c*o = -(1 + 0.325) -p;f. (6.12) 



The remaining parameters of the polynomial part start at 0(j) ) and we get them from (5.5) 



( p!25| ) and ( ji^ ). Their numerical values and the related subthreshold parameters Cij in natural 



units (chosen in such a way to make the comparison with the polynomial parameters easy, i.e. we 
take cio and coi in units of M'^/F'^ and C20 in units of F""*, cf. ( 4.16| )) are shown in Tables [| and 
I for the two sets of 0{p'^) LEC's. 

All the considered parameters are strongly sensitive to the variations of the LEC's. For instance 
the parameter (3 varies with the LiS as 

A/3 = (17.OAL1 - 2.8AL3 + 9.lAi4)M2. (6.13) 

For the LEC's Q we get (3 = 0.90 x lO^^M^. Using the set [|5) we get a drastic change 

A/3 = -1.58 X lO-^M^. (6.14) 
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Let us turn to the "doubly dangerous" observables represented by the scattering lengths now. For 
the s-wave we obtain from ( 4.15| ) and the LEC's Q 



ao = ^^^^(1 + 0.683 + 0.378- 0.625 + 0.527) 



(1 + 0.963) = 11.0 X 10^^ (6.15) 



247rF2 Mr, + M, 

Here the individual terms in the first line represent the current algebra result, the correction stem- 
ming from the 0{p'^) contributions to the parameters a and ui, the next-to- leading correction to 
induced by the expansion of the denominator and the correction induced by the dispersive part of 
the amplitude 0(s), in this order. This result confirms the expectations about a bad convergence of 
the chiral expansion for the observables which are connected to the threshold values of the ampli- 
tude - even if the polynomial NLO corrections were small, which are not, the dispersive part would 
be still as large as 50% of the leading order term. 

The sensitivity to the O(p^) LEC's is illustrated in Table ||. The p-wave scattering length then 



starts at 0{p'^), we get the values in the last column of the table from ( 4.15 ). 

When comparing our results for the scattering lengths with those of , quoted in the Intro- 
duction, we can see a large discrepancy 

Uq = 1 .2 X 10 Ui = —5.2 X 10 . (6. lb) 

The difference is produced by a different set of 0(p'*) LEC's, the alternative treatment of the 
in the denominator and by another form of the unitarity corrections - the authors do not use a 
matching with a dispersive representation. Taken these distinctions into account, we get more 
consistent numbers (with our inputs for the masses and decay constants): uq — 7.0 x 10~^ and 
ai — —5.0 X 10~^. As we can see, a slightly different treatment of the standard chiral expansion 
may lead to a significant shift in results. 

6.2 Resummation of vacuum fluctuations - basic properties 

In the resummed case, the free parameters are X, Z, and r together with the remaining LEC's 
Li, L2 and L3 and the direct and indirect remainders da---S^ and Sp^^ SppMp- Because is 
experimentally not known with enough accuracy, we also have to fix how to treat the observable 
^GMO which was introduced to eliminate the LEC L7 using the bare expansion for F^M'^. Let us 
remind that our definition of Aqmo follows the ref. [|| , where it is based on the good observables 
FpMp instead of Mp and differs from that originally defined in [^. One possibility is to treat 
Agmo a-s an additional independent parameter. The another, similarly to the treatment of -Fi, 



in the denominators of (4.14), (4.15) and (4.16) discussed earlier, is to use a (resummed) chiral 
expansion of F^ inserted to Aqmo for the numerical estimates, i.e. to insert the following exact 



algebraic identity into ( 5.13 ) (cf. Appendix g for details) 



FlMl^GMo = FlMl - AMlFl + M^il - Sp^^ (^4F|.(1 ~ 5f^) - F^{1 - SpJ 

-Ml (^1) (j;^(0) - 2(r + 1) J],^(0) + (2r + 1) j;,(0))) . (6.17) 

This generates the indirect remainder Sp^ in a nonlinear way. 

Before doing a more detailed analysis, let us first illustrate the numerical sensitivity connected 
with the subtleties of the definition of the bare expansion. As we have discussed in Section 3, there 



is still some freedom how to define the amplitudes entering the dispersive part of the G^,, (cf. (3.3C 



3.31 ) and ( 3.32 , 3.33 ) and also how to treat the masses inside the chiral logs. Based on general 
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considerations it was argued Q , that in the latter case the different prescriptions should not make 
much difference. Nevertheless, it might be interesting to test this assumption numerically in our 
concrete case and also to check what is the numerical influence of the varying amplitude definition. 

In Figure |^ we plot the comparison of 



10 



2 47tF„^ (M^+M„) /M„^ ao at Z = 0.9,r = 20 



various definitions of the dispersive part of 
the amplitude using the scattering length oq 
as an example, i.e we illustrate its sensitivity 
on the various versions of the unitarity cor- 
rections. The cusps on the full line, which 
uses the strict chiral expansion with the un- 
physical choice of the 0{p^) masses in all 
Jlj, originate in the conflict of the physi- 
cal masses used for the on-shell outer legs 
and the unphysical location of the thresh- 
olds. This illustrates the fact that the orig- 
inal strict chiral expansion is unsuitable for 
realistic physical predictions and its redefi- 
nition into a bare one is necessary. The dot- 
ted line shows the "minimal" physical mod- 
ification of the strict expansion by means of 
insertion physical masses into all J^j . While 
The horizontal line shows our standard NLO prediction. the "minim al" version and the unitary choice 

(3.32, 3.33) give numerically almost the same 
result, the difference between the these two and the third possibility (3.3C, 3.31) is up to ^ 0.3 a^^. 
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Figure 1: Comparison of the numerical impact of the 
various forms of the dispersive part on the scattering 
length ao . The full line represents the strict chiral expan- 
sion, dotted, dashed and dash-dotted lines the "minimal" 
modification, (3.3C, Wsl) and {3.3i, 3.3t) respectively. 



at Z = 0.8,r = 20 
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Figure 2: In this figure we illustrate the sensitivity of the "good" variables a and (5 to the treating of 
the chiral logs. The full line corresponds to the 0{p^) masses in all ,J[j{0)'s, while doted and dashed lines 
to the physical masses either in all Jlj (0) 's or only in Jlj (0) 's originating from the unitarity corrections. 
Horizontal lines are standard NLO predictions. 

Figure || shows the dependence of the polynomial parameters a and (3 otlY — X/Z for Z = 0.8 
and r = 20, using the the various possible treatments of chiral logarithms in the bare expansion. 
The results demonstrate that the difference might be numerically important in some range of Y . 
For a the various possibilities do not differ drastically in comparison with the value of a itself, on 
the other hand the differences become comparable with a^^ at y ~ 0.5. As we have discussed 
in Section 2, in the region of small Y the case with 0{p^) masses in the tadpoles only should not 
differ drastically from the case when all the masses are physical. However, the convergence to the 
common value at y = is rather slow and in the intermediate region of Y the difference of this two 
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cases for a is ~ 0.5 oP^ in a relatively wide interval. Keeping 0{p^) masses also in the unitarity 
corrections produce instabilities for y — > as expected. The parameter /3 (which starts at 0{p'^)) 
is even much more sensitive. 

In the following numerical analysis we take a pragmatic point of view and fix the bare expansion 
in such a way that the comparison of the resummed and standard reparametrizations remains as 
simple as possible, i.e. we insert physical masses into Jlj{0) and define the amplitude according to 
(|I28| , and 



6.3 Numerical comparison of the resummed and standard reparametrization 

Within the standard we have an 0{p'^) prediction for X, Z, r and Aqmo based on the standard 
formal 0{p'^) chiral expansion ( ]a1| ) and ( [A^ ), cf. Appendix Using the LEC's from @ and 



||4q| , we get numerically the following central values, which should confirm the self-consistency of 
the standard chiral expansion scheme 
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As we can see, while expectations are fulfilled in the first case, there is a considerable shift when using 
the O(p^) fitted constants. These numbers, moreover, should be taken with some caution, because 
they originate in the expansions of the "dangerous " observables and can be therefore plagued with 
large 0{p^) remainders as well as with strong sensitivity to the 0{p'^) LEC's^"'^. In the above table 
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Table 3: The values of the the polynomial parameters a and /3 and the related subthreshold and threshold 
parameters near the standard reference point. For Aqmo we take either the standard value (Aq^iq ) or the 
resummed prediction described (Aqmo ) in the main text. 



Let us now illustrate the relationship of the resummed and standard approach using the ob- 



servables from Subsection 3.1 



std 



For the "good" observables a and /? we can expect that the numerical values of X 
and A^ffQ, with Li, L2 and L3 taken from B3] for definitencss, should produce numbers consistent 



with the first row of Table |l| when inserted into (5.13, 5.14). The results for the various possibilities 
how to approach the standard predictions for a and /3 (which is independent on Ag]\jo) within 



^'^As it was analyzed in detail in g], the actual values of X"^'^ and Z"^'^ are strongly sensitive to the values of the 
LEC's Lg and L4 connected with the vacuum fluctuation of the ss pairs, the same is true for the sensitivity of r**'' 
and Aq'^jq to Lg and L7. This causes large error bars to be attached to these values. Nevertheless, in the following 
we take these central values as a reference point for an illustrative numerical comparison of the two versions of the 
chiral expansion. 

^■^Though the difference between the values of r''*'* and r*''*'* is within the standardly expected accuracy of the 
0{p'^) approximation, note, however, that for r*^*'* the O(p^) correction is much larger than in the first alternative 
: 39.4). 
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Figure 3: The dependence of the parameters a and j3 onr, X and Z is plotted, one of the parameters 
fixed at its standard reference value in each figure. The dashed horizontal line shows the standard values 
from the first row of the Table ^ the full circle depicts the corresponding resummed value at the standard 
reference point [r^*'', X^*'', Z"*"^] . The error bars represent the 10% uncertainties from the direct and indirect 
remainders added in quadrature. In the first row, r is fixed at r"*'*, the filled are highlight the dependence 
on Z between Z = Z"*** (solid line) and Z = 0.5 (dotted one). Similarly, in the second row X = X**'', the 
filled area shows the dependence on Z again. Z is fixed at Z"^''' m the last row, the solid line shows the case 
with r = r"'"* while the dotted the one with a lower value r = 15. 



the resummed version of are summarized^^ in the Table The last row corresponds to the 
resummed treating of Aqmo explained above. The dependence of the central values of a and P 
on the parameters r, X and Z in the broader vicinity of their standard values is illustrated in 
Fig.^. These results can be interpreted as a consistency of both variants of reparametrization for 
good observables near the standard reference point X'**'', Z''*'* and r^***, where the predictions of 
the resummed version almost coincide with the standard results^"*. This coincidence together with 
the working hypothesis about the controllable remainders of good observables within the resummed 

■^^In this and the following tables in this subsection we ignore the uncertainty stemming from the remainders and 
Li, i = 1,2,3 and give only the central values (assuming the central values of the remainders to be zero). 

■^"^As a rule, the point Jf"*'*, Z"*'^ and r"*'^ cannot give the best coincidence with the standard values in all cases, 



- 31 - 



F^Vm„2 lO^Cio at r=r 























T 






1 




=1 





















0.2 0.4 0.6 0.8 

X 

Coo /Coo"'*' at X=X='=d 



16 18 20 22 24 26 28 30 
r 

Coo /Coo'''*' at Z = Z=td 




0.2 0.4 0.6 OA 
X 




0.2 0.4 0.6 O.f 

X 

F^Vm^2 lO^Cio at X=X=^="^ 



16 18 20 22 24 26 28 30 

r 

F^Vm„2 lO^Cio at Z = Z=^='* 























1 


H 


] 1 


1 T 

























.2 0.4 0.6 OA 
X 



Figure 4: Dependence of the subthreshold parameters cqo and cio related to the polynomial parameters a 
and p. The figures are in one-to one correspondence to those in Figure H. 



reparametrization scheme confirms again a self-consistency of the standard expansion based on 
the assumption X ~ 1, Z ~ 1 and r ~ r2. Away from the standard reference point, however, 
the standard reparametrization might be dangerous in the sense that the difference between the 
standard and the resummed prediction diverges rapidly and the importance of the standard 0{p^) 
remainders might therefore increase considerably. 

For the "dangerous " observables like Cij we cannot a priori expect coincidence of both expansion 
even near the standard values of X, and r due to the different treatment of the denominators, 
which contain large Oip^) corrections and are not expanded in the resummed case. Comparison of 
both approaches is illustrated in Table ^ (with the same treatments of A^mo ^ above). Table ^ 
and in Fig. ^. For the dispersive part we use the prescription ( 3.32| , 3.33| ), which differs from the 



corresponding standard contributions of the unitarity corrections to Cy for X = Z = 1 and r — r2 
by a factor F^/F^ w 0.6. This is reflected by the values of those Cy that start at O(p^) (cf. Table 
H, ^). Namely in this case the contribution of the polynomial part is reduced near the reference 
point roughly by the same factor with respect to the standard value (which includes only the first 
term of the expansion of the denominator). On the other hand, coo is compatible with the standard 



the reason can be understood e.g. by having a closer look on the resummed reparametrization of /3 (cf. (5.14)). In 
order to reproduce the dependence of /3 on L4 satisfactorily, we need Z = and r = r''", on the other hand to 

reproduce the chiral logs we need rather X/Z = 1 and r = r2. This can explain why /3 approaches the standard 
value best for X = Z = Z^". 
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Figure 5: Dependence of the scattering lengths ao and ai related to the polynomial parameters a and fi. 
The figures are in one-to one correspondence to those in Figure [J. 



value, because of the large OijP') contribution, tiny dispersive contribution and the fact that within 
the standard reparametrization of the bare expansion also the second term from the expansion of 
the denominator is taken into account. 

Let us now proceed to the "doubly dangerous" 
obscrvablcs oo and oi. These are related to the values 
of the amplitude at the threshold and receive there- 
fore large contribution from the dispersive part of the 
amplitude. While ao is reproduced well at the stan- 
dard reference point, ai (which starts at the NLO) is 
off the standard value roughly by a factor 0.6 from the 
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Table 4: The values of the subthreshold pa- 
rameters C20 and coi related to the polynomial 
parameters 7 and uj at the standard reference 
point. 



same reasons as for the Cy parameters. The depen- 
dence of this observables on X, Z and r is depicted 
in Fig. |5[ 



6.4 The role of the remainders 



Up to now we have not discussed the uncertainties of the observables calculated within the resummed 
scheme. They are connected with the direct and indirect remainders as well as with the LEC's L^, 
i = 1, 2, 3. As the first illustration, we have added the error bars stemming from the remainders to 
the central values of the various observables depicted in the figures |^||. These illustrate the rough 
estimate of the remainders 6 ^ (30%)^ ^ 0.1 as suggested in |^ and adding the uncertainties in 
quadrature. 

In more detail, at the standard reference point X = X^*'^, Z = Z'^^'^ and r — r^*'^, using ( 5.15 



5J6| ) and (U, |]T|), we numerically get for the corresponding variations (to the first order in the 



remainders) 
Aa 

+ (3.38ALi + O.56AL3) X 10^ (6.18) 
A/3 = [(0.69(5/3 + 2.MSf^ ~ 20.52(5fJ + (IT.OALi - 2.8AL3) x 10^] x lO^^Af^. (6.19) 

This reveals strong sensitivity on both (5s and LEC's. Assuming again the typical size of the 
remainders to be (5 ~ 0.1 and adding all the uncertainties in quadrature (for AL^ we take the error 
bars form E3]) we obtain rough (over)estimates 



= v/1.932 + 1.192 = 2.27 (6.20) 



|A/3| = ^2.062 + 5.962 ^ 10-3^2 = 6.31 x 10-='m2, (6.21) 

where the first number under the square root represents the contribution of the remainders while 
the second accumulates the uncertainty from £1.3. Though these numbers are a little bit more 
optimistic than those in the last row of Table |l| (note that the latter originated purely from the 
uncertainties of AL^ and did not include any estimates of the higher order corrections to a and 
/3), it is clear that, without more restrictive information on the remainders (and i^, i = 1,2,3)2^, 
the predictive power of xPT is reduced considerably in the case of ■kt] scattering even for "good" 
observables. In other words, small remainders are not a guarantee of an equivalently small final 
uncertainty. In what follows, we therefore try to gain some additional information outside the 
(resummed) xPT expansion to get further estimates of the size of the remainders. 

The sources of the remainders are twofold: on one hand there are the unknown terms of the 
pure derivative expansion, on the other hand the contributions coming from the expansion in the 
quark masses. We try to get estimates for both of them from different sources, namely using the 
resonance estimate for the first type as well as independent information from generalized X-PT for 
the second. 



6.5 Resonances estimate of the direct remainders 

In order to partially estimate the derivative part of the higher order corrections to the chiral 
expansion, we use the assumption, that the process under consideration is saturated by the exchange 
of the lowest laying resonances, the interactions of which can be described by the Lagrangian of the 
resonance chiral theory (i?xT). The leading order Lagrangian of RyT was originally formulated in 
the seminal paper Q and appHed to 7r?7 scattering in js^. To this process, only scalar resonances 

^^As already (discussed, the explicit dependence on these constants could be eliminated by means of reparametriza- 
tion similar to those for L^, i = 4, . . . , 6 using further experimental input e.g. form ife4 decay. The price to pay is 
to introduce additional remainders connected with observables used for such a reparametrization. 
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as well as rjg, — 770 mixing contribute. Our result for the amplitude agree with js^ (cf. Appendix 
0), which we can rewrite in terms of the resonance contribution G^,^ to G-^jj in the form 

Gf^is, t; u) = af + f3ft + + Lof{s - uf + AG^,(s, t- u). (6.22) 

The polynomial part with the coefficients (in what follows, M5 and Ms^ are the octet and singlet 
scalar mass respectively, Cd, Cm, Cd, c„i and dm are the couplings defined in |^6| ) 



2 ( ^dCm 'iCdCjn \ , ^ f 



"'[fMi-iiC)- '"lit I"' ' " 

m!!-/: - J,) - m! ^ + 4 m: + ) (6,23) 



.(4) _ SCdCd , 4 C2 / c^c„ C<iC„\ / o 2 o 2\ 

(4) _ _ ^ 25) 

4'' = 7^ (6.26) 
^ 3M| ^ ' 

gathers the complete O(p^) resonance contribution (here we can recognize the resonance saturation 



of the leg's in (4.23)-(4.26)). This part of the amplitude is already included in our resummed 
version of xPT, either explicitly through the LEC's Li. . . L3 or implicitly using the reparametriza- 
tion in terms of the masses, decay constants and parameters r, X and Z. On the other hand, 
AG7r,7,i{(s, t; u) can be formally understood as an infinite sum of the higher order corrections in the 
(purely) derivative expansion, summed up to 



AG^^^fl(s,i;-u) = 



(jvf 2 , ^) [^dit - 2Ml) + 2?™ Af, j [cd{t ~ 2M2) + 2Zm M,j 



2 



2 u 



2 



2 t 



3M|(M|-t) 



o 2 

Cdit-2M^) + 2c„,M^ 



X [ crf(i - 2Af2) + 2c™(2 M,', - mI) 



Of course, this does not exhaust all of the possible higher order corrections (note e.g. that the 
resonance Lagrangian we use contains only the leading order interaction terms with one resonance 
field and chiral building blocks of the order 0(p2)), nevertheless we can use it at least as a rough 
estimate of the effect of higher orders of the derivative expansion. This is in some sense a procedure 
opposite to the usual resonance saturation; instead of LEC's we "saturate" the remainders by means 
of sewing together the resummed chiral expansion G^^^^ (without remainders) with the resonance 
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chiral theory writing the fuU RxT amphtude as 

G^^^is, t- u) = G^7(s, t- u) + AG« (s, t- u). (6.28) 

and identifying G^^^ with the fuU xPT amphtude, the remainder being AG^^. Under this as- 
sumption, we can derive the fohowing higher order contributions to the direct remainders from 
AG^^ 

16/ o'^\ ( o2 

A/3ii = /S-^^ = UdMl - Cm MA ( CrfM^ - Cm M„ 

8 1/ , o2\ / o2 o2 



(6.29) 



, o4 o2 o4 



3M|(M|-I]^^) 3M|(M|-I]^^) 3M|(M|-E^^)2 

and similarly for 6!^ and 5^- (see Appendix 0) . Note that the dependence on X and Z is exclusively 
through the ratio Y = X/Z here. 

One can notice that there are two distinct features of this procedure as compared to the usual 
LEG saturation. First, there is no need to fix a saturation scale, which is the result of "saturating" 
the renormalization scale independent remainder instead of the scale dependent LEC's. And second, 
as the resonance contribution are resummed to all chiral orders, the resonance poles are explicitly 
present in our result, as can be seen in ( 6.27 ), ( 6.29| ) and ( 6.3C ) as well as the formulae for 5^ and 



5^ in the Appendix |j[ 

For the rough numerical estimates we use M5, M^^, M^^ and the couplings Cd, Cm, 'cd, 'cm, dm 
from [Q. This gives at the standard reference point X = X**'', Z = Z'**'' and r = r^*** 

5f 1.00 (6.31) 
I35f = -0.15 X IQ-^M'^^, (6.32) 

which represents roughly 55% and 20% correction to values in the first row of the Table ^ respec- 
tively. The dependence of 5^ and 5^ OT^Y = X/Z and r is depicted in Figure |[ The effect of the 
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resonance remainder estimate on the parameters a and /? in a wider range of the X, r and Z is 
illustrated in the first column of Figure analogous plots for 7 and uj are in Figure ^ As can be 
seen, these results suggests the conclusion that the derivative part of the expansion could in some 
cases produce higher order remainders with much bigger value than 10%. 
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Figure 7: Polynomial parameters 7 and u depending on Y . Horizontal dashed line: standard 0{p'^) and 
central RxPT value. The result with resonance remainder estimates is shown by the sohd line. 



6.6 Generalized xPT 

In the previous subsection we have tried to estimate the contributions to the remainders generated 
by the derivative expansion. The resulting expressions ( |6.29 ) and (6.30) could, however, gather 
only terms of at most the second order^^ in the quark mass expansion due to the lowest order 
resonance Lagrangian used. Also, the indirect remainders have not been included in this way as 
there is no contribution to them in this simplest approach. For the appraisal of the importance 
of the missing terms we therefore need additional information. One possibility might be to use 
a resonance Lagrangian with additional terms of higher chiral order suited for saturation of the 
0{p^) leg's and/or to go to the next-to-next-to leading order in the chiral expansion, this is, 
however, beyond the scope of our paper. 

Instead we try to get some flavour of the size of the effect by means of comparison of our previous 
results with generalized xPT, which was originally designed to handle the badly convergent quark 
mass expansion in the case X <^ 1 and therefore also includes terms which correspond to higher 
orders in the standard chiral power counting. 

In Subsection 5A , we have already rewritten the generalized expansion of the parameters a and 
/3 (as well as that of the masses and decay constants in the Appendices ^ || ) in the "resummed" form 
( |5l3Hn6| ) by means of splitting the "standard" remainders into the "nonstandard" extra terms 
^ioops^^-j g^^^ ^'^^(/i) originating in GxPT and the unknown part 6'^^^^. Therefore, neglecting 
the latter, the sum (5'°°P*(/i) + 6^'^{ii) could be in a sense interpreted as the rough estimate of 
the contribution to the standard remainders stemming from the higher orders of the quark mass 
expansion. 

While (5'°°P'*(/i) are known, S'~^'^{^) depend on the unknown LEC's of the GxPT Lagrangian (cf. 
Appendix ^). We therefore set S'^^{iS) = at a fixed scale /i and by varying this scale in 5^°°'p^(ii) 
from jjL = to /i = Mp we can get some information on the contribution of the unknown LEC's 
(note that 5^°°p^{ii) +5'~^^{ijl) is renormalization scale independent). We apply this procedure both 
to the direct and indirect remainders. 

The usual way of handling the generalized xPP expansion is to neglect the unknown remainders 
SGxPT^ We can repeat the considerations from the previous subsection and partially appreciate 

■^^Note that the physical masses in (6.2E), (6.3C) originate in the derivative expansion. 
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Figure 8: Polynomial parameters a and 13 depending on X and Z for traditional and low value ofr. The 
dotted line shows the central value for Z = 0.9, the dashed one is the same for Z = 0.5. The error bars 
correspond to the 10% estimates of the remainders. Left column: resonance estimate, filled areas highlight 
the 0{p^) and higher corrections to the amplitude generated by resonances (lighter for Z = 0.9, darker for 
Z = 0.5). Right column: results with combined resonance and GxPT estimate of remainders. Filled areas 
show the scale dependence ~ M,, — Mp), lighter for Z — 0.9 and darker for Z = 0.5. 



them using the resonance estimate. In order to avoid double counting, we have to further modify 
the resonance contribution to the remainder ( 6.27 ) subtracting terms of the order O(p^) within the 
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generalized power counting in the same way as it was done in the previous subsection (c.f. ( |6.22| )) 



This combined GxPT and resonance estimate of the remainders is ihustrated in Figure the 
right column shows the result in the case of the polynomial parameters a and f3. The effect of 
the unknown GxPT LEC's is estimated by their scale dependence. The lines closer to the central 
i?xPT results with neglected remainders are the ones at the scale /i = Mp , i.e. the constants are 
set to zero at the usually chosen scale. The filled grey areas then show the change when the LEC's 
are set to the difference when moving from the scale ^ = Mp to M^. Admittedly, this assigns quite 
arbitrary numbers to the LEC's, so the uncertainty should be viewed as a rough estimate which 
can go both ways. The result can be interpreted as being quite consistent with the 10% estimate 
of the remainders, though clearly exceeding it for some range of the free parameters X, Z and r. 

As for the parameters 7 and u, because their contribution in the polynomial expansion is 
quadratic in the Mandelstam variables, the GxPT estimate docs not contribute here. 

7. Summary and conclusions 

In this paper we have studied the properties of various variants of the chiral expansion on the 
concrete example of the nrj scattering with special attention to the possible reparametrization in 
terms of the physical observables. We have tried to illustrate several issues in detail, namely: 

• the necessity of carefully choosing a class of "good" observables for which the condition of 
global convergence is believed to be satisfied in the sense that the 0{p^) and higher remainders 
are small and under control. 

• the necessity to carefully define the bare expansion of "good" observables. Here we have 
concentrated on the requirements dictated by the exact renormalization scale independence 
as well as the exact perturbative unitarity. As we have shown, both these requirements can 
be met by means of sewing together the strict chiral expansion in terms of the LEC's with the 
dispersive representation for the corresponding Green function. Nevertheless, the resulting 
bare expansion is not yet defined uniquely; one has to fix the way how to treat the chiral 
logs and also the 0{p^) amplitudes entering the dispersive integrals. Though the difference is 
formally of the same order as the remainder itself, we have found that, it might be numerically 
significant in some region of the free parameters. 

• the properties of the standard chiral expansion, based on the potentially "dangerous" reparametriza- 
tion of the bare expansion implicitly assuming X, Z 1 and r r2. In this case we have 
established strong sensitivity of the observables for Tryy scattering on the O(p^) LEC's; this 
plagues the standard prediction with large uncertainty. In the case of L4 and ig this also 
means strong sensitivity to the vacuum fluctuations of the Is pairs and therefore to the devi- 
ation from the standard scenario with X, Z ^ I. The unusually large absolute values of the 
NLO corrections as well as large variations achieved for most of the observables (including 

the "good" ones) when moving from the 0{p'^) fit of the LEC's Li Q| to the 0{p^) based 
fit |45) might be interpreted as a signal of the importance of the NNLO corrections within 
the standard chiral expansion. This seems to be also supported by the sensitivity of the 
NLO contributions to their form when expressed in terms of the physical masses and decay 
constants {i.e. how the 0{p^) relations like e.g. the Gell-Mann-Okubo formula are used). 



^G'^Sis^t^u) = AG 




— 4 d (M, + M2)- Mr, t) 




(6.33) 
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• properties of the "safe" reparametrization and resummation of the vacuum fluctuations. We 
have confirmed that, for the "good" obscrvables, the rcsummcd and standard values coincide 
near the standard reference point [X Z**'*, r^*'']. Under our working hypothesis, which 
assumes the "good" observables to be accompanied with smaU and controUable remainders, 
this can be interpreted as a consistence of the standard O(p^) chiral expansion of "good" 
observables in the sense that the potentially large higher order remainders are in fact small. 
On the other hand, in most cases of the "dangerous" observables the standard and resummed 
values do not meet at [X'**'^, Z**'*, r (typically for those of them that start at 0(p'*)). 
Though this might indicate that the standard expansion is convergent less satisfactorily in 
this case and the higher order remainders might be important here, the difference between 
the standard and resummed values lies within the estimated uncertainty of the resummed 
prediction. Away from the standard reference point, however, we have established that the 
central values diverge substantially from those of the standard approach even for the "good" 
observables. This is a signal that, unless X, Z ^ 1, the higher order remainders of the 
standard chiral expansion might be huge in comparison with LO+NLO value. Though this 
feature does not exclude the possibility that in this case the standard remainders might be 
saturated by the NNLO corrections, it could be nevertheless interpreted as an indication of 
the instability of the standard chiral expansion. 

• the role of the remainders within the resummed approach. We have found a strong sensitivity 
of the observables connected to ttt] scattering to the higher order remainders. This might 
reduce the predictive power of this approach, unless additional information on the actual 
size of the remainders is available. We have tried to make an independent estimate of the 
remainders using the simplest version of the resonance chiral Lagrangian as well as making a 
comparison with GxPT. Both these estimates seem to be in accord with the rough expectation 
6 ~ 10% for the remainders only in some range of the parameters. For some observables 
and some corners of the parameter space, they can be substantially larger. Of course, the 
convergence properties of the bare expansion deserves further investigation by means of going 
to the NNLO, which is, however, beyond the scope of our article. 

As a final remark, we would like to comment on the outlook of future usability and practicality of the 
outlined procedures in a broader picture. One can easily argue, that even if the 0{p^) standard xPT 
was considered inadequate, proceeding one step further in the chiral expansion, as was reviewed 
e.g. in |Q, is the next logical step to do. Large 0{p^) remainders could alter the precision of the 
results, but hardly their overall validity. On the other hand pushing the resummation an extra 
order brings difficulties with the exact algebraic reparametrization. 

However, at standard xPT 0(p^) order one encounters the ultimate practical limit of effec- 
tive theories, the growing number of LEC's. Additional information beyond core xPT has to be 
gathered, usually resonance saturation is used. In this sense, the resummed xPT framework in 
the present form makes a barrier between the explicitly manageable part of the chiral expansion 
and the remainders which need to be estimated at the right point. If e.g. exactly the same RxT 
Lagrangian to 0{p^) was used, the method of incorporating resonances outlined above would in 
fact retain more information about the complete RxT theory than usual Oljfi) LEG saturation. 
Also, explicitly calculable contributions like double logs could be included into remainders similarly 
to one-loop generalized xPT terms in this paper. Thus, the possible advantages of resumming the 
vacuum fluctuations could be exploited with only moderate additional computational effort. 
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no. LC 527) and by the EU Contract No. MRTN-CT-2006-035482, "FLAVIAnet". 
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A. Standard chiral expansion of the parameters X, Z and r 

Here we summarize the formulae leading to the standard values of X, Z, Y and r used in Subsection 
6.3. Using the standard reparametrization rules explained in Subsection 5.1, we get up to the NLO 
order in terms of the 0{p'^) LEC's 



X 



std 



^Std 



l-^(32(Lg(r2 + 2)+LD 

+3J;,(0) + (r2 + l)JU(O) + i(2r2 + l)j;,(0) + ^\l^J'^ 
1 - ^ U {Ll{r, + 2) + LI) + {V2 + 1) J]^k(0) + 4J;„(0) - 



r2 



167r2 



std 



r2 



2F1 



8(21; - L,){r2 - 1) - \{2r2 + 1) j;,(0) + JIM %^ 



For r we can also use an alternative expression based on the chiral expansion of 



(A.l) 



2FI 

+ i(2r2 + 1) j;,(0) + + 1) J],k(0) - | j;,(o) + 



Agj\/o has the following standard chiral expansion 
Ml{r2 - 1) 



967r2 



^GMO 



2F2 



(32(2L'7 + Lg)(r2-1) 



1 



+ -(2r2 + 1) j;,(0) + (r2 + 1) J],^(0) - 3j;,(0) + 



5(r2 - i; 
487r2 



(A.2) 



(A.3) 



B. Chiral expansion of the r) decay constant 

For the bare expansion of the "good" observables Fp we rewrite the standard formulae in the form 



Ft - F' 



-^0 



- 



1 + 



Bom 

^0 



16i^(M)(r + 2) + 16L^(m) + {r + 1) J],^(0) + 4j;,(0) + ^(r + 5) 

1 



(B.l) 



3 

+ 2 



16L^(M)(r + 2) + 8L^(/.)(r + 1) + -(r + 1) J]^^(0) + 



+ -(2r + l) ( j;,(0) + 



167r2 



1 + 



Bofh 



16L^(M)(r + 2) + ^L^(M)(2r + 1) + 3(r + 1) ( J],^(0) + 



167r2 



(B.2) 



(B.3) 



Within the standard xFF, the 0{p^) parameters Bq, Fq and r are expressed using inverted expan- 
sions of the observables Fp, Mp, as explained in subsection 5.1. This yields the standard formula 



for F2 



F' 



F^ 



Ml /16 



p2 



^■LlinXr^ - 1) + (r2 + 1)J],^(0) - 2j;,(0) + ^(r^ - 1) 



(B.4) 
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with a potentially large remainder 6p . Numerically, with Ll{Mp) taken from |Q we get 



(B.5) 



On the other hand, the "safe" reparametrization in terms of r, X and Z gives 

(j;,(0) - 2(r + 1) J],^(0) + (2r + 1) JIJO)) 



f2 _ ui2 



(B.6) 



which is valid as an exact algebraic identity^''. 

Following the GxPT procedure outlined in Section 5.4, after identifying the corresponding 
leg's in both approaches 



B^fh 

^0 



(B.7) 



and defining the remainders using the physical masses inside the chiral logs, we can use the exact 
formula (B.3) and write the remainder 6f^ within GxPT as 

Here 

F^5%°^{^^) = 3m^ (r + 1) (^o (r + 1) + 2(r + 2)Z^) (^r^^iO) + ^) (B.8) 

is the extra loop contribution and the contribution of the counterterms from the O(p^) 

GxPT Lagrangian renormalized at scale /x 



^{2Ai +A2 + 4A3 + 2Bi - 2B2)(1 + 2r2) 



(B.9) 



+3(^4 +2Bi){l + -r^) -4Cf(r-l)2 + 2D^(r + 2) (2r+l) 

^GxPT ^ remainder, which is exactly independent on the renormalization scale. 

Analogously, for the GxPT formula for we have, besides the substitution (B.7) to (B.l), to 
insert 

F^Sf^ = FXrif^) + Fo%'! ""^(f^) + FXf^, (B.IO) 



^'^This identity can be also rewritten as 

4Fi(l - 5^^) - f2(1 ^SfJ- 3F2(1 - SpJ = J) (j;,(0) - 2(r + l)jp,-^-(0) + (2r + l)j;„(0)) . 

Within the standard approach, the parameters on the r.h.s. of this identity can be expressed to the order 0(p*) 
in terms of the physical observables and it is interpreted as a O(p^) sum rule 

4F|. - F2 - 3F2 = m2 ( J;^(0) - 2(r2 + 1) J]f^(0) + (2r2 + 1)J;^(0)) . 



This gives 
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where the loop and counterterm contribution are now 

FXTit^) = 8™2(Ao + {r + 2)Zl)) (^j;,(0) + ^) 



+fh' (r + 1) (Ao (r + 1) + 2(r + 2)^5' ) J]f;^(0) + 



167r2 



Ai + + 2^3 + (A4 + 2B4)(1 + ir^) 



+Bi - S2 + 2D^ [r + 2)] . 
FinaUy, we have the expression for where the remainder is replaced with 

and the loops and counterterms contribute as 

1 



(B.ll) 



(B.12) 



FlS^'it^) = 3^2(^0 + (r + 2)^0^)) j;,(0) + 



167r^ 



+ |m2 + 1) (Ao (r + 1) + 2(r + 2)Zo^) I J],^(0) + 



167r2 



+m2(Ao (2r2 + l) + 2(r - l)^^^^ + (r + 2)(2r + 1)Z^) (^j;,(0) + ^) 



(Ai + Bi)(r2 + 1) + {A2 - 2B2)r + 2(^4 + 2S4)(1 + -r^) 



+2i:>^ (r + 2) (r + 1)] . (B.13) 

In order to reparametrize the GxPT bare expansion in terms of the masses and decay constants, 
we can proceed as follows. Because the exact identity (B.6) is valid independently of the version of 
xPT, we can also use it in the generalized case, provided we rewrite the remainders according to 
(B.9, B.IO) and (B.12). This step eliminates the LEC's ^ and ^. Collecting the chiral logs together 
we have 



p2 _ pi2 

7/ TT 



1 



1 



167r^ 



167r2 
1 



F^^/4f^{^^) (B.14) 



where the 0{p^) masses are given by ( ^.19 ) and 



(B.15) 

The last step consists of replacing the LEC's Fq, Aq, Zq and Zq with the first term of their expan- 
sion in terms of the masses and decay constants as described in Subsection ^.4| . This corresponds 
to a further redefinitions of the generalized remainders. 



C. Dispersion representation of the ttt] amplitude 

For the dispersive representation of the amplitude we need the S— and T— channel discontinuities 
at O(p^). In the following subsections we give a list of the relevant 0{p^) amplitudes G^^-''^'^*-' and 
Qi^)ij^^f ai4(j O(p^) discontinuities discGQ-' corresponding to the different intermediate states ij. 
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c.l 5— channel discontinuities at O(p^) 
• nr] intermediate state 



discGo^Cs) = 2 ^ ^ — - j ^ (C.l) 



• KK intermediate state 



1 ^ o2 o2 

1 ^ o2 o2 o2 



V327ry 

1 o2 o2 o"^ ^F^ 

--{2Mk-M^-M,)?^ (C.2) 



C.2 T— channel discontinuities at O(p^) 
• TTTT intermediate state 



discGr''=°(s) = 2a(s) (^)' i mI [{s - \mI) + \ Mj^ (C.3) 



■qr] intermediate state 



KK intermediate state^^ 



4 uiv 3 3 

2 o 2 o 2 

+ -(Mk + M,)] 



g(2)xoxO(^+^-)^^^_,^o ^ ±iF2[(3s-2M|.-2M2) 



^*Let us note 

G'=°[s,t;u) = ^<5<'''G°''(s,t;u) = -s/ZG{s,t;u) 

v3 
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o2 2 o2 

+(2M^--Mk)] 



2 o 2 o 2 

+ 3(M^ + M,)] 

X [(35 - 2Mi - 2M2) + (2 M?,)]] ^ (C.5) 

D. The scalair bubble 

In this appendix we summarize the formulae for the scalar bubble, defined as 



{2-kY (fc2 - Ml + iO)((fc - g)2 - M2 + iO) 

„2\ 



= -2A^ + rpQ{q'). (D.l) 
Here, as usual 



(D.2) 



and Jpq{s) = Jpq{0) + Jpq{s), where 



1 M|,ln(M|/M^)-M^ln(M^/M2) 

and Jpq{s), sometimes called Chew-Mandelstam function, can be expressed by means of once 
subtracted dispersion relation as 



^ dx \^'\x,Ml,Ml) 1 

'{Mp+MqY 



167r2 ./rA/f„4-»^„^2 x x x-s 



The explicit form of Jpq{s) reads 

7pg(.) = f 2 + ^ m ^ - ^ m ^ + 2(^-(^^-^Qn,,^(,) 1, ^£4f)_l\ ^ 

"^^^ ^ 327r2 \^ s Ml ApQ s ^"^^ ' apQ{s) + 1J' 

(D.5) 

where 

ApQ = Mj,-M^ 

2 I n#2 



SpQ = M| + M, 



. - (i\/p + 4M,,J/g 



In the limit Mp — > Mg we get 



1 / /Vf2 
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E. O(p^) constants L^-Lg in terms of masses and decay constants 

In this Appendix we summarize the formulae used in the text for the reparametrization of bare 
expansions of "good" observables. We use the abbreviated notation (5. £- 5.1^ ). From the bare 
expansion of "good" variables and F|. we obtain 



2' " "r + 2 

Ml X 



4(r + 2)(r- 1) Z 

(r + 2)(r-l) 



[(4r + l)j;,(0) + (r-2)(r + l)JU(0) 



(2r+l)j;,(0) + 

lU7f- 

(E.l) 



2(r + 2)(r- 1) 



. 2 



4 Llif,) = ^F^ryW 



Af2 X '\(r — 1 

+ T7-^^[5JU0)-ir + l)JKK{0)-i'2r + l)JU0)- ^ 



4(r-l)Z^ ""'^ ' y-KKv-; V j-^vvy-' iQ^2 



(r-1) ■ 

In the same way, from the expansion of FpMp we get 



o 4 1 F'^M'^ 

4 r + 2 



(E.2) 



72(r - l)(r + 2) ' ^''^'^^'^ + + '^^'^ " '^^^-^'^ 



+ (2r + l)(r - 4) j;,(0) + ig^l'^^^ ] (^.3) 



f>/^,5F,M. [(r + 1)^] - 4F|-Mi(5f^,,M.. 
4(r2 - l)(r + 2) 



I 4 



4m! = ]F^A4e{r) 



3(1 + r)F^M^SF,M, + (2r2 + r - 1)F^M^Sf^m^ ~ 8rF^Mj,SF^M^ 

8(r-l)2(r + l) 

F^W^rfr) 

4 



(E.4) 



^ ^ [9j;,(o)-3(r + i)j;,^(o)-(2r + i)j; ro)- ^ 



(E.5) 



24(r-l)V^/ ^ '^'^'^ ^ J-KKK-j y j-^n^y-j -^g^2 

2FIM15f^m^ - (r + l)F,^M2jj.^M. 
2(r2 - 1) 

F. Lagrangian of GxPT to O(p^) 

Here we give the traditional form of the GxPT Lagrangian. In the following formulae 

X = M + s + ip 
VU = dU - i{v + a)U + iU{v - a) 

X = dx~'-i{v + a)x + ix(w - a). (F.l) 



-46 - 



Up to the order O(p^), the Lagrangian can be spht into the 0{p^), 0{p^) and O(p^) parts 

£ = £2 + £3 + £4, (F.2) 

where 

£„= £(''3'^) (F.3) 

i+j+k—n 

and k) indicates the number of derivatives, x sources and powers of Bq respectively. Then for 
O(p^) we get 

£(2,0,0) ^ El{^^u+WU) 

£(°'^'°) = ^(Ao((c/+xf + (x+t/)^) 

+ {U+x + X^Uf + XiU+X - X^Uf) ■ (F.4) 

At the order O(p^) one has 

+P3((x+t/f -(c/+xf)(x+t/-c/+x) 
+p4((x+c/)' + (t^+x)')(x+t^ + c^+x 
+Mx-^u + c/+x)(x+x) + P6(x+t^ - u+x?{x^u + C/+X) 
+P7(x+c/ + c/+x)') 
£(2,0,1) ^ ^^^w^v^cz+v^c/) 

£(0,2,1) ^ :^(^(l):4^^(f;+^)2^(^+f^)2 

+ 5(i)4(c/+x + x+c/)' + s^'K{u^x - x^uf) 

£(0,1,2) ^ ^B^,5(^'Hu+x + X^U). (F.5) 
For the O(p^) Lagrangian, the building blocks are 



£(4,0,0) ^ Ll{V^U+V^'U)^ + L2{Vf,U+V^U}{V''U+VU) 
+i3(V^t/+V^C/V,C/+V"C/) 

(U+F^^UF^,) + H, {F^,F^^^Ftf,F^"f') 
£(2,1,1) ^ El^(^5Wliy^u+W^U{x^U + U+x))+5^^%^^.U+V^'U){x+U+U+x)) 

£(2,0,2) ^ £;M42)^v^[^+V''i7) 

£(2,2,0) ^ El{A^l^y^u+V^U{x+x+U+XX^U)) 
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+A2{{\/^U+)Ux+{^''U)U+x) 

+A4(V^C/+V^C/)(X+X) 
+Bi{V^,U+WU{x^Ux+U + U+xU+x)) 
+S2(V^C/+XV''C/+X + X+V'f/x+V^t/) 
+B4(V^[/+V''[/)(x+?7x+f/ + U+xU+x) 
+t?f (V^f/X+ + xV^f/+)(V''C/x+ + XV''C/+) 

+c|(v^x+c/ + t^+v^x)(V'^x+f/ + t^+Vx) 

+C^3^(V^X+C/ + U+W^x){^^U+x + X+V^!7) 

+cr(v^c/x+ - xv^c/+)(v^t/x+ - xv^;7+) 
+c^l'(Vpx+(^ - (:^+v^x)(v^x+c/ - c/+v^x> 
+C^3''(VmX+C/ - t/+v^x)(v^i/+x - x+v^c/) 
+D^(V^t/+V'^C/(x+;7 + U+x)){x^U + U+x) 
+D^(V^C/+V^C/(X+J7 - U+x)){x^U - u+x)} 

+H2 (v^xV^x^) 

£(0,4,0) ^:^{^^^(^+^)4^(f^+^)4^ 

+£^2(x+x(x+(:^x+c/ + c/+xc^+x)) 

+E3{x+xU+XX^U) 

+Ff{x+Ux+U + U+xU+xf 
+Fi{{x+Uf + iU+xf){x+U + u+x) 

+if (x+x(x+f/ + u+x)){x^u + u+x) 

+Fi{{x+U)' + {U+x)'){x^x) 

+Ff{x+Ux+U-U+xU+x)^ 
+Fi'{ix.+U)' + iU+x)'){x+U + U+X) 

+Fi'{x'^xix^u~u+x)){x^u-u+x) 

+Jf^(x+x)(x+c/ + (:^+x)' 
+i^/^((x+c/)' + iu+x)^){x+u - u+x) 

+Fr{x+x){x-'U-U+x)' 
+Ff''{{x+Uf - {U+x)') 
x{x+U-U+x){x^U + U+x) 
+H3 (XX+XX+) + Hi (XX+)' 

£(0,3,1) ^ ^^s^i)j,^^{x+U)' + {U+x)')+S^'^p,{{x^U + U+x)x^x) 

+^'''P3((x+c/)' - (c/+x)')(x+t^ - u+x) 

+5(i)p4((x+c/)' + (c/+x)')(x+t^ + u+x 

+s^'%{x^u + c/+x)(x+x) + s^'^Mx^u - u+x)'{x^u + u+x) 

+5(^)-p,{x+U + U+x)^) 
£(0,2,2) ^Ei^(^5(^)^^i^^u+x)'^{x-'Uf) 
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(F.6) 



In fact, identifying Fq with the Goldstone boson decay constant and Bq — Yj / where S = —{uu)q 
(in the chiral hmit), we have 



5« = = 0. 

As usual, we can also resume the powers of Bq already at the Lagrangian level and write 



i+j—n 



with 



and denote 



k) 



= + BoS^^^Ao + BlS^^^Ao + 



Z, 



S,P 



+ B,5^^^ZI-^ + Bl5^^)X-'' + 



^S.P 



■S,P 



p,, = ft + Bo^'i^p, 



(F.7) 
(F.8) 

(F.9) 



(F.IO) 



etc.. These LEC's without the bars are then used in the main text. Note that while the O(p^) 



parameters Zq and the 0{p ) LEC's ^ ,^ are renormalization scale independent, the renor- 



malized resumed parameters Zq'^'^, Aq and run with /i in the same way as 16{Bo / Fq)'^ L^^g^ 

and 8Bq/FqL4^5 within the standard xPT. 



G. CoefRcients of the dispersive part of GxPT amplitude 

In these formulae as well as in the following two appendices, the masses Mp are the generalized 
O(p^) masses given by ( 5.19 ). 



a M2 



Mi = 



a^^(4M2 - Ml) 



2[mBo + Bm^A^ + 2m'^Zl{^r + 4) - 8m2Z^(r - 1)] 
Am\T^ ^1){A^ + 2ZP) 

2mBQ + Wui^Aq + 4m'^Zf^{r + 8) 

2 16 

-mBo(l + 8r) + —m^Ao{l + 8r^) 

o o 

32 

i + 41r + 32r^) + —m'^Z^{r ~ l)(4r - 1) 
3 



-m'Z^i 



6{Aa + 2Z^)m^{r + l) 
2 

-[mBa{l + 3r) + m^Aa{3 + lOr + lOr^) 

O 

+2m^Z^{6 + 19r + Ur^) + 16m^Z^r{r - 1 



(G.l) 



H. Parameters a — uj within the generahzed xPT 

Here we summarize the formulae in terms of the decomposition of the remainders. For the parameter 
a we write 

'5"=^i°°^ + 3^e(A^)+e^^. (H.l) 
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For the counterterm contribution we get 

= [81 pi+ P2 + (80 - 64 r - 16 r^) ps 

+ (100 + 64 r + 34 r^) p4 + (2 + r^) ps 

+ (96 - 96 r) pe + (l44 + 288 r + 108 r^) pr] 



1 



-(i3i - B2)S^r, + 2i?^M2(r - 1) - 2Cf M2(r - 1) - -i^^^[S,^(5r + 4)] 
- 2B4[3M2 + M2(2r2 + 1)]] 

+ im2 [256£i + 16^2 + (256 - 256r2) + if (32 + ler^) 

o 

+ Ff (256 + 32Gr2) + ^5^^(192 - 32Gr + leOr^ - 32r3) 
+ F2^(240 - 216r - 24r^) + Fff^(32 - 32r + 16r^ - 16r^) 
+ ^3^(16 - 8r - 8r^) + Fi{16 + lOr + lOr'') 

+ Fg^^(32 + 40r + ler^ + 2Gr3) + F/^(384 - 160r - 2567-^ + 32r3) 

+ ^2^(400 + 234r + 74r^) + F^^{576 + 720r + 480r2 + 168r3)] (H.2) 



and the loops contribute as 



^F^M^J^-f = ^ I [m' (3Bom + 64^om2 + 2Z^m\l5r + 32) - SZ^rh^Sr - 8))] 



-6i?o'm2}(j;,(0) + ^ 
+ ^ |[M^ (Bom + 2Aofn^{r + 8) + 2Z^m^{15r + 8) - 8Z^fn^{3r - 2))] 
-B^fh'ir + l)} (^j;,^(0) + ^) 

i I [M^{Bofh + 32Aom^ - 16Z^fn^{5r - 2) + 2Z^m^{4:lr + 16)] 



+ 1 |[m' + Am\3Ao - 4(r - 1)Z^ + 2{2r + l)Z^)f - 4Bo'm^} J;,(0) 



+ 1 {[^M^ - ^(r - l)fh\Ao + 2Z^)f - ^Bo^m^l J^^(O) 
+ 1 |[m' + Afh\3A^ - 4(r - l)Zo^ + 2(2r + l)Zl)] 



X [-2Ml + -M2 + 10m2(Ao + 2Zl)\ - 2Bom{3Bom - 2M^) } j;^(0) 



+ 1 [{mI + im\3Ao - 4(r - 1)Z^ + 2(2r + 1)Z^)] 



X - -M2 + m\{8r^ + l)Ao + 8r(r - 1)ZP + 2{2r + ifZ^) 



1 

vi;- 

■\Blffi\8r + l)\r^^{Q) 



-i [[-2Ml + 2Ml + 8(r + l)fh\Ao + 2^0^))] 
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X [-6M,; + 6M;; - + -{r + l)m^3rAa + 2{r - 1)Z^ + 2{2r + 1)Z^] 

-2{2Bom - M2)(lBom(4r + 1) - 6M^)^ J^kkW- (H.3) 

In the same way we have for (3 

(iSp = + m'FM^i^,) + /3<5«>^^^, (H.4) 



where 



Vif') = gK^f + D'){2r + 1) + 2B,{r' + 2)] (H.5) 
+ 1 [Wl + 4a^(3Ao - 4(r - l)Z(f + 2(2r + - 2Bom| JIM 



+ _(^ + l)m2(3Ao(r + 3) + ^Zl{r + 5) + 2(r - 1)ZP)] 

-[^i?om(2r + 5) - GM^ - GM^lj J],^(0). (H.6) 
For the remaining two parameters the corresponding decomposition of the remainders 

7<50(m) = T^-r'^^^lM) + ^?^'J^o + 7'^?''''^ (H.7) 

c.(5^(m) = c.<5^'°°^^(m) + m'n'C^(M) + ^•5^'^''^ (H.8) 

is trivial, i.e. 

= S^^ii,) = = <55^(a*) - 0. (H.9) 



I. Generalized xi'T contributions to the bare expansion remainders for 
the masses 



The expressions for ^, ^ can be obtained from the exact algebraic identities (E.2) after identification 
(5.24) and using the representation (B.IO) and (B.ll) for the remainder of and ( B.12| ) and (B.13) 



for the remainder of i^^. In the same spirit, Aq, Zq and Z^ can be expressed using identities (E.5) 



and the following remainders 



F^^M^Sf^^m,, = F^^M^S'^Ji, (^) + Fo^m^^gfM, (a^) + ^'M,^^^,^^;^, (I.l) 
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where 



5fImM 



then 



"FkMk (m) 



Ml {Worn + IQAofn? + 2Z^w?{ir + 16)) - 6S^m^ 
1 

167r2 



X JIM + 



+2 



M|: (Bom + 2Aofh^{r + 2) + 2Z^m^{2,r + 4)) - B^TO2(r + 1) 
1 



X ( J^k(O) 

2 



1 

+ 3 



167r2 



M^{Bom + SAofh"- - SZ^m^r - 1) + 2Z§m\5r + 4) 

1 



~B2m'(2r + l) 



./,%(0) + ■ 



(1.2) 



m[9pi + /)2 + 2p4(10 + 4r + r^) 
+P5(2 + r^) + 12p7(4 + 4r + r^)] 

+2m^[8£i + 2£;2 + 8Ff (2 + r^) + Fi (9r + + 20) +Fi {4: + r + r^) 
+2Fi (2 + r^) + 4^5^^ (r + 2) (r^ + 2r + 6) + 2i?^g^^ (r + 2) (2 + r^)], (1.3) 



(Bom + Aofh^ir + 5) + 2Zf m2(r + 6)) - 2B^m' 
1 



X j;.(o) 



167r2 



+ 2 



(Bom + 3^om^(r + 1) + 2Z^m^{3r + 4)) - B^fh^{r + 1) 



X J^x(O) 



1 



167r2 



1 

+ 12 



(5BoTO + Aom2(17r + 5) + 8Z^m^{r - 1) + 2Z^m2(13r + 14)) 
1 



-^Bo^m^(2r + l) 



JL{0) + 



167r2 



(r + 1) 



(1.4) 



-m[3pi(l + r)(l + r + r^) + p2{l + r^) + 6p4 (r + 1) (2 + 2r + r^) 

+P5 (r + 1) (2 + r^) + 12p7 (r + 1) (r + 2)^] 
+m2[2Si (1 + rf (1 + r^) + -Ea (1 + rf (l - r + r^) 

4^3(r^-l)^ 

+4i^f (1 + rf (2 + r2) + F^^ (1 + r) (8 + 9r + 9r2 + Ar^) 
-Fi (1 + r) (4 - r + + 2r'^) + Pf (1 + rf (2 + r^) 
+4i^5^'^ (1 + r) (2 + r) (4 + 3r + 2r2) 
+2^^6^'5 (1 + r) (2 + r) (2 + r2)] 



(1.5) 



and 



Fv^v^fJi, (/«) = (^om + 8 ^0 - 8 Z^ (r - 1) + 2 m^ (4 + 5 r) ) - 2B^fh'' 

1 



X ^;.(o) 



167r2 
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+ 1 [m|- (Bom (l + 4r) + 2Aom2 (2 + r + 8r2) +8fh^Z^ {r ~ I) (2r-l) 
+2 a^Zf (4+ 15r + Sr^)) - Sgm^ (4r + 1) (r + 1)] 



167r2 



1 



M^{Bom (l + 8r)+8Aom2 (l + Sr^) +16m2z^ (r - 1) (4r-l) 



+2 a^Z^ (8 + 41r + 32r2)) - -B^m^ (8 r + 1) (2r + 1) 



X j;,(o) + 



167r2 



*F;m, (m) = 3^19^1 (1 + 2r') + p2 (1 + 2r3) + 16p3 (r - 1)' (1 + r) 
+2p4 (10 + 8r + ITr^ + lOr^) 

+P5 (1 + 2r) (2 + r^) + 16pe (2 - 3r + r'^) + 12p7 (2 + r)^ (1 + 2r)] 

+^m2[8Bi (1 + 2r4) + 2^2 (l + 2r4) 

+8Ff (r^ + 2) (2r2 + l) + WPf (r^ - l)^ 
+Fi (20 + 13r + 37r3 + 20r'*) + 12^2^ {r^ + r + l) (r - 1)^ 
+^3^ (4 + 5r + Sr^ + 4r*) + AF^ (r^ + r + l) (r - 1)^ 
+2Fi (r^ + 2) (2r2 + l) 

+ l2Fi"^ (r + 2) {2r^ + 3r^ + 2r + 2) + 8Fi^ (r^ + 2) (r - 1)^ 
+2i^g^^ (2r + 1) (r + 2) (r^ + 2) + 4^6^^ (r^ + 2) (r - 1)^ 
+ 16i^/^ (r + 2) (r + l)(r-l)']. 



(1.6) 



(1.7) 



J. Resonance amplitude and remainders estimates 

Here we give the contribution to the amphtude p8[ | related to the resonance exchange, derived 
from the leading order Lagrangian of RxT (we have confirmed this expression by independent 
calculation) 

Gb{s, t- u) = ^J^p—^ {cd{t - 2Ml) + 2Zm {cait - 2M^) + 2Z„, 

Cp o / o ^ o2\ g2 Q 2 / Q 2 Q 2 

i^M^ + M„j + 4^ i^M^ - 

1 / o 2^ 2 

Cd(s - Af2 - M^) + 2Crn 



V 



2 1 



3Mj-u 



3M| 



(^Cd{u - Ml - Af2) + 2Cra M^ 

Cd{t - 2Ml) + 2c™ mI^ (^Cd{t ~ 2M^) + 2c„(2 m] - Mj^ 



(P o / o o2 

J)! 



-16^ ( - Mj . (J.l) 
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The resonance estimate of the remainders 5?- and 5,^ are 



7 

= (^CdM^ - Cm M^j i^CdM^ - Cmi2 Mr, - M ^) 



-I nf -I iJl yr 



f ^ _ o 2 o 2 \ 

CdS^^ - Cm{M^ + M^) (J.2) 



'Si 

■31 



_ 1 C^S,^ 4 I 4 CjCm o ^ , . 

iMl{Ml-^^r,)^{Ml-T.^r,f^^{Ml-T.^r,?^' ^ ' 
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